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Abstract 

We study Compton scattering in systems with A=l and 2 using chiral pertur- 
bation theory up to fourth order. For the proton we fit the two undetermined 
parameters in the 0(Q 4 ) 7p amplitude of McGovern to experimental data in the 
region u, y|ij < 180 MeV, obtaining a x 2 /d.o.f. of 133/113. This yields a model- 
independent extraction of proton polarizabilities based solely on low-energy data: 
a p = (12.1 ± 1.1 (stat.))!^ (theory) and j3 p = (3.4 ± 1.1 (stat.))!^ (theory), 
both in units of 10~ 4 fm 3 . We also compute Compton scattering on deuterium 
to 0(Q 4 ). The 7d amplitude is a sum of one- and two-nucleon mechanisms, and 
contains two undetermined parameters, which are related to the isoscalar nucleon 
polarizabilities. We fit data points from three recent 7d scattering experiments with 
a x 2 /d.o.f. = 26.6/20, and find a N = (13.0 ± 1.9 (stat.))^ (theory) and a P N that 
is consistent with zero within sizeable error bars. 



PACS nos.: 13.60.Fz, 12.39.Fe, 25.20.-x, 12.39.Pn, 21.45.+V 



1 Introduction 



Compton scattering on the proton, neutron, and deuteron provides a window on the inter- 
nal dynamics of these strongly-interacting systems. As this is an electromagnetic process, 
the photon-target interaction can be treated using well-controlled approximations. At the 
same time Compton scattering is governed by a different set of target structure functions 
than, say, electron scattering, and so probes complementary aspects of the strong nuclear 
force. For example, in the case of the proton, quark models indicate that even at quite 
low energies 7p scattering can reveal distinctive information about the charge and current 
distributions produced by the nucleon's quark substructure (UI21- 

The deuteron represents a different theoretical challenge. Its structure is governed by 
the nucleon-nucleon (NN) interaction. A number of models of this interaction exist which 
very accurately describe experimental data on iViV scattering at laboratory energies be- 
low 350 MeV. However, these successful approaches differ significantly in their underlying 
assumptions. They range from models which are almost purely phenomenological [3] to 
multi-boson-exchange models |H E] that make detailed assumptions about how QCD 
plays out in the NN system, e. g. the role of the pomeron E|. Clearly, NN scat- 
tering alone cannot completely constrain the form of strong interactions at low energies. 
Deuteron Compton scattering is a way of probing the deuteron bound-state dynamics 
generated by different models. Furthermore, any theoretical treatment of this reaction 
must incorporate not just photon interactions with the target via the one-body Compton 
process but two-body mechanisms as well: the input to a calculation of the 7d amplitude 
includes the jN amplitude plus Compton-scattering two-body currents. These two-body 
currents should be constructed in a way that is consistent with both the single-nucleon 
Compton dynamics and the NN potential used to generate the nuclear bound state. 

Over the last several years a framework [7j based on effective field theory (EFT) [S] has 
been developed that allows such consistency. Importantly, EFTs also satisfy the symmetry 
constraints of QCD while making only minimal theoretical assumptions. This framework 
can be applied to the Compton processes we are interested in, provided that the photon 
energy is well below the chiral symmetry breaking scale, A x ~ 47rf n ~ M ~ m p . In this 
regime photon-nucleon and photon-deuteron scattering are governed by the approximate 
chiral SU(2)i x SU(2)r symmetry of QCD. The spontaneous breaking of this symmetry 
to isospin SU(2) V results in the existence of three Goldstone bosons, which are identified 
with the pions (7r + , 7r~, and tt°). 

An EFT can be built from the most general Lagrangian involving pions and nucleons 1 
which is constrained only by approximate chiral symmetry and the relevant space-time 
symmetries. In this EFT, known as chiral perturbation theory (%PT), pions interact 
through vertices with a countable number of derivatives and/or insertions of the quark 

1 The simplest form of the EFT results if the Delta isobar is integrated out of the theory, although 
this tends to limit the range of energies over which data can be described. 
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mass matrix. Therefore S'-matrix elements can be expressed as simultaneous expansions 
in powers of momenta and pion masses over the characteristic scale of physics that is not 
included explicitly in the EFT, A x . A generic amplitude can be written as: 



where Q represents the external momenta and/or the pion mass, T v is a calculable function 
with ^v(l) of order one, C is an overall normalization factor, and v is a counting index. 

%PT has proven remarkably successful, with a number of single-nucleon processes 
including 7p scattering [TO] ITT| IT2"] — now computed to several nontrivial orders. In xPT, 
as in any EFT, the price to be paid for the absence of model assumptions is a lack of 
detailed information about physics at distance scales much shorter than the wavelength 
of the probe. This price translates into the presence in xPT amplitudes of certain unde- 
termined interaction strengths, or "low-energy constants" (LECs). These constants can 
be estimated using models of the short- distance physics, but in the purest form of %PT 
they should be fitted to experimental data. 

There has been intense recent effort dedicated to extending %PT to processes with 
more than a single nucleon p]. The fundamental difficulty in making this extension is 
the treatment of the effects of nuclear binding. At momenta comparable to the pion mass 
the few-nucleon EFT first proposed by Weinberg still provides a paradigm for the 
application of EFT to nuclear systems. Use of this few-nucleon EFT has developed to 
the point where computations of processes involving two nucleons with any number of 
pions and photons can be carried out at levels of precision comparable to those achieved 
using ;\TT in the single-nucleon sector jHj. The present paper is part of this ongoing 
program [13 EE! • 

During roughly the same period over which EFT has developed as a tool for analyzing 
low-energy QCD dynamics, experimental facilities with tagged photons have made pos- 
sible a new generation of Compton scattering experiments which probe the low-energy 
structure of nucleons and nuclei. An extensive database now exists for Compton scat- 
tering on the proton at photon energies below 200 MeV with this database being 
considerably enhanced by the recent data taken in the TAPS setup at MAMI by Olmos 
de Leon et al. [IEj. In the nuclear case the past decade has seen a set of experiments that 
establish — for the first time — a modern data set for Compton scattering on the deuteron. 
Data now exist for coherent 7c? — > 7c? from 49 to 95 MeV ^Hl 123 l2*Tj . and for quasi- free 
7 d -> 7pn from 200 to 400 MeV [221 123 121 ■ 

These experimental and theoretical developments come together in studies of Comp- 
ton scattering on the A = 1 and A = 2 system in xPT. Nucleon Compton scattering 
has been studied in xPT to 0(Q 3 ) HO] and 0(Q 4 ) [UJ H2J- To 0(Q 3 ), there are no 
undetermined parameters in the 7p amplitude, so xPT makes predictions for 7p scatter- 
ing. These agree reasonably well with low-energy differential cross-section data [3 EE] ■ 




(1) 
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At the next order in the chiral expansion, 0(Q 4 ), there are four undetermined param- 
eters, two for 7p scattering and two for scattering. These counterterms account for 
short-range contributions to nucleon structure. Minimally one needs four pieces of exper- 
imental data to determine their values, and once they are fixed xPT makes improved, 
model-independent predictions for Compton scattering on protons and neutrons ^2J- Be- 
low we use low-energy proton data in order to obtain the two short- distance parameters 
associated with the proton, and so make a completely model-independent determination 
of the xPT 7p amplitude up to fourth order in small quantities. 

In the A = 2 system, the amplitude for coherent Compton scattering on the deuteron 
was computed to 0(Q 3 ) in Ref. There are no free parameters at this order. The 

corresponding cross section is in good agreement with the Illinois data ^H] at 49 and 
69 MeV, but underpredicts the SAL data [H] at 95 MeV. The calculation of Ref. [T5] 
also yields cross sections which agree well with the more recent Lund data [21]. In the 
present paper we extend this calculation to 0{Q A ). At this order coherent 7d scattering 
is sensitive to two isoscalar combinations of the four free parameters in the 0(Q 4 ) 'jN 
amplitude. We will show that fitting these two parameters to experiment yields 7d cross 
sections which are in good agreement with the existing data for momentum transfers and 
photon energies below 160 MeV. A summary of earlier 0(Q A ) ;\RT results for both 7p 
and 7d scattering was presented in Ref. ^B]- The current paper contains a considerably 
more detailed discussion. Because of errors in the fits to 7d data reported in Ref. [IB] the 
central values found here for the isoscalar polarizabilities differ significantly from those 
given in Ref. ^B] (see Section IB~4l for a full explanation). 

Very low-energy nucleon Compton scattering (with uj significantly below m^) has been 
a focus of particular theoretical interest. At these energies, the xPT 'yN amplitude can 
be expanded in powers of uu/m w . This yields a T-matrix that is a simple power series in 
u, with coefficients which are functions of %PT parameters. In the nucleon rest frame, 



Here the ellipses represent higher powers of energy and momentum as well as relativistic 
corrections, and e (e') is the polarization vector of the initial- (final-) state photon and 
k (k') is its three-momentum. The first term in this series is a consequence of gauge 
invariance, and is the Thomson limit for Compton scattering on a target of mass M and 
charge Z\e\. Meanwhile, the coefficients of the second and third terms are the nucleon 
electric and magnetic polarizabilities, a and (3. To 0(Q 3 ) in %PT they are completely 
given by pion-loop effects [TU] : 
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with qa the axial coupling of the nucleon and f w the pion decay constant. We emphasize 
that the polarizabilities are predictions of x?T at this order: they diverge in the chiral 
limit because they arise from pion-loop effects. In less precise language, ;^PT tells us that 
polarizabilities are dominated by the dynamics of the long-range pion cloud surrounding 
the nucleon, rather than by short-range dynamics, and thus should provide a sensitive 
test of chiral dynamics. 

At 0(Q 4 ) the four undetermined parameters alluded to above contribute to the polar- 
izabilities [TT|. They can be understood as representing short-distance (r <C effects 
in a and (3. The results from our fits to proton Compton data provide the first completely 
model-independent determination of proton polarizabilities. These results, together with 
an analysis of various remaining sources of theoretical uncertainty, are given in Section 0] 
below. 

While data from Compton scattering on hydrogen targets can be used to extract the 
proton polarizabilities a p and f3 p , the absence of dense, stable, free neutron targets requires 
that the neutron polarizabilities a n and f3 n be extracted from scattering on deuterium (or 
some other nuclear target). Coherent Compton scattering on deuterium depends on the 
polarizability combinations = (a p + a n )/2 and (3n = (/3 p + /3 n )/2 through interference 
between the polarizability pieces of the 7iV amplitude and the proton Thomson term. 
Thus ot j\f and (3^ can be extracted from nuclear data, but to do so requires a consistent 
theoretical framework that cleanly separates single-nucleon properties from multi-nucleon 
effects. In the long-wavelength limit pertinent to polarizabilities EFT provides a model- 
independent way to do exactly this [7J |HJ E] , and thus we would argue that (or at least 
some systematic EFT) is an essential tool in the quest to obtain and f3^ from deuteron 
Compton data. Our x?T analysis of 7d data in Section El facilitates a model-independent 
extraction of a at and (3n, as well as a quantitative assessment of the uncertainty associated 
with the freedom to choose different NN potentials when generating the deuteron bound 
state. 

Our results for both 7d cross sections and extracted values of a at and [3n agree qualita- 
tively with potential- model treatments of 7d scattering at these energies [2H1 I27J 12H] • The 
different potential-model calculations all yield similar results if similar input is supplied — 
in particular in the form of comprehensive meson-exchange currents. Typically though, 
the one- and two-nucleon Compton scattering mechanisms in potential models are not 
derived in a mutually-consistent fashion. The calculations of Refs. |2ol 1271 125] also tend 
to incorporate significantly more model assumptions about short-distance physics than 
are necessary in an EFT approach. 

We note that calculations of Compton scattering on the deuteron in EFTs other than 
the one employed here exist in the literature. For example, in Ref. [29J, an EFT where 
pions are integrated out is employed. This EFT is under good theoretical control and is a 
very precise computational tool, but it is limited in its utility to momenta well below the 
pion mass. In Ref. j3U] a different power counting |3T] from ours was used in which pion 
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exchange is treated in perturbation theory. While this power counting is well-adapted to 
processes at momenta below the pion mass, it is now known to break down at quite low 
energies 

The rest of the paper is structured as follows. In Sectional we summarize those aspects 
of xPT which are directly relevant to our calculation, including the power-counting scheme 
and the pertinent terms in the chiral Lagrangian. In Section El we explain how this theory 
is applied to Compton scattering on the proton up to 0(Q 4 ), reviewing the results of 
Ref. ^2]- In Section0]we fit the two unknown parameters in the 0(Q A ) x?T 7p amplitude 
to 7p scattering data with oj, \/\t\ < 180 MeV. In Section we turn to the A=2 system, 
computing the Feynman amplitudes that contribute to photon scattering on the two- 
nucleon system at 0(Q 4 ). In Section |H1 we present results for differential cross sections 
for Compton scattering on unpolarized deuteron targets. We discuss higher-order effects 
and the extent to which a determination of neutron polarizabilities from elastic 7d is 
possible. Finally, in Section Owe summarize our work, and identify future directions for 
the theoretical study of this reaction. 

2 Baryon Chiral Perturbation Theory 

In this section we briefly discuss the effective chiral Lagrangian underlying our calculations 
and the corresponding power counting. 

2.1 Effective Lagrangian 

QCD has an approximate chiral SU(2)l x SU(2)f> symmetry that is spontaneously broken 
at a scale A x ~ 1 GeV. Below this scale there are no obvious small coupling constants in 
which to expand. Effective field theory is the technique by which a hierarchy of scales is 
developed into a perturbative expansion of physical observables. Here we are interested 
in processes where the typical momenta of all external particles is p <C A x , so we identify 
our expansion parameter as p/A x . 

In a system with broken symmetries this technique is especially powerful. When a 
continuous symmetry is spontaneously broken there are always massless Goldstone modes 
which dominate the low-energy dynamics. If the chiral symmetry of QCD were exact, 
p/A x would be the only expansion parameter. However, in QCD SU{2)l x SU(2)r is 
softly broken by the small quark masses. This explicit breaking implies that the pion has 
a small mass in the low-energy theory. Since m n /A x is then also a small parameter, we 
have a dual expansion in p/A x and m 7T /A x . 

We will limit ourselves to the region where ~ p < — M. (M and Ma stand for 
the nucleon and Delta-isobar mass, respectively.) In this regime we can keep only pions 
and nucleons as active degrees of freedom, and the low-energy EFT is without an 
explicit Delta-isobar field. We take Q to represent either a small momentum or a pion 
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mass, and seek an expansion in powers of Q/A x . We implicitly include the mass scale 
of the degrees of freedom that have been integrated out — Ma — M, m p , etc. — in A x . 
One might reasonably expect that the convergence of the theory will improve upon the 
inclusion of an explicit Delta field [UnilS]- Work on nucleon Compton scattering along 
these lines is ongoing (351 EE] • 

Although chiral symmetry is no longer a symmetry of the vacuum, it remains a sym- 
metry of the Lagrangian. One can show |S7] that a choice of pion fields is possible, in 
which all pion interactions are either derivative or proportional to m\. The symmetries 
of QCD still permit an infinite number of such interactions though, and so it is necessary 
to order the strong interactions according to the so-called index of the interaction [T3] . 
For an interaction labeled i this is defined by: 

A, = di + fi/2 - 2 (4) 

with di the sum of the number of derivatives or powers of m n and fi the number of 
nucleon fields. Since pions couple to each other and to nucleons through either derivative 
interactions or quark masses, there is a lower bound on the index of the interaction: 
Aj > 0. We write the effective Lagrangian as 

oo 

£ = ££ (A) . (5) 

A=0 

So far this only accounts for strong interactions. The electromagnetic field can enter 
£, either through minimal substitution or by the addition of terms involving the electro- 
magnetic field strength tensor. The former simply has the effect of replacing a derivative 
by a factor of the charge e. And in either case we now have another expansion: one in 
the small electromagnetic coupling a em = e 2 /47r. However, since we work only to leading 
order in this expansion it is convenient to enlarge the definition of di so that it includes 
powers of e, and then continue to classify interactions according to the index Aj defined 
by Eq. ®. 

The technology that goes into building an effective Lagrangian and extracting the 
Feynman rules is standard by now and presented in great detail in Ref. [§]. Here we list 
only the elements necessary for our calculations. 

The pion triplet is contained in a matrix field 

where f„ is the pion decay constant, for which we adopt the value f„ = 92.4 MeV. Under 
SU(2) L x SU(2) R , S transforms as £ -> LY,B) and f as f L£W = U£I&\ here L(R) 
is an element of SU{2) L (SU(2) R ), and U is defined implicitly. It is convenient to assign 
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the nucleoli doublet N the transformation property iV — > UN. With Z = (1 + t%)/2 we 
can write the pion covariant derivative as 



D tl Z = d tl Z-ieA t t[Z,Z}. (7) 



Out of £ one can construct 



V, = ^(d^-ieA^ + ^-ieA^) (8) 

A, = ^{dp-ieA^Z-tidr-ieAftZ)?], (9) 

which transform as V„ -> UV^W + Ud^W and A„ -> fL4 M f/t under 5C/(2) X x 5C/(2) K . 
is used to build covariant derivatives of the nucleon, 

D^N = (fy + V„)N. (10) 

(Ref. j2Hl uses where we have used and for our A^.) The electromagnetic field 
enters not only through minimal coupling in the covariant derivatives, but also through 

f, u = e(t}Zi + ZZt})F^ (11) 

where F^ u = d^A u — dyA^ is the usual electromagnetic field strength tensor. 

Because the nucleon mass is large, M 3> Q, it plays no dynamical role: nucleons are 
non-relativistic objects in the processes we are interested in. The field N can be treated 
as a heavy field of velocity v in which on-mass-shell propagation through exp(iMv ■ x) 
has been factored out. In the rest frame of the nucleon, the velocity vector t> M = (1,0). 
The spin operator is denoted by 5 M , and in the nucleon rest frame = (1/2)(0,<t). 
This formulation of is sometimes called heavy-baryon chiral perturbation theory 
(HB X PT). 

With these ingredients, it is straightforward to construct the effective Lagrangian. As 
will become clear later, in order to calculate Compton scattering to 0{Q i ) we need to 
consider interactions with Aj < 3. 

The leading-order Lagrangian is 

= ^ 2 Trp M Et^E) + i/>^Tr(S + Et) 

+iN\v ■ D)N + 2g A N\A ■ S)N 
C 

~Y [3N^N N^N + N^SN ■ N^SN] + ... (12) 

with g A = 1.267 the axial-vector coupling of the nucleon and Cq a parameter determined 
from NN scattering in the 3 Si channel. The next-to-leading-order (NLO) terms are 
contained in 

= jM Ni {- D2 + (v D ) 2 + 2 *9A{v-A,S-D} 



-'-[S^S" 



+2M Cl m2Tr(E + E f ) + (8Mc 2 - g 2 A ) (v ■ A) 2 + 8Mc 3 A 2 } N + ..., (13) 
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where k v = k p — n n and k s = k p + n n are parameters related to the anomalous magnetic 
moments of the proton, k p = 1.79, and neutron, K n = —1.91; and C\, c 2 , and c 3 are 
parameters that can be determined from ttN scattering. We use c\ = —0.81 GeV -1 , 
c 2 = 2.5 GeV -1 , and c 3 = —3.8 GeV -1 jSH] (compare with similar values obtained from 
NN scattering [40J). The ellipses in Eq. (J13|) indicate that we have not written down 
pieces of which are not relevant to our computation of 7iV and 7d scattering. 
At sub-sub-leading order (N 2 LO): 

£( 2 ) = 7Ti^ up(T F F 

32n 2 U 3 1X1/9(7 

-y [SN*NN* (-D 2 + (v ■ D) 2 ) N + N ] SN ■ N ] S (-D 2 + (v ■ D) 2 ) N + h.c] 

-y [N\S ■ D)NN\S ■ D)N + h.c] + ... (14) 

where the first term is from the chiral anomaly, with £0123 = 1, and the other two terms 
are corrections to the two-nucleon potential in the deuteron channel, with the coefficients 
G% and C e determined by fits to iViV scattering data. 

In Eq. (I14j) we do not explicitly display terms which are needed for our computation 
below and have coefficients that are determined by Lorentz invariance. These coefficients 
are dimensionless numbers times 1/M 2 , and the relevant numbers can be extracted from 
Ref. [2E|- A fixed-coefficient 7^ seagull from which enters the NLO 7iV amplitude 
stems from the 12th and 13th terms of Eq. (3.8) in Ref. |38j . while the '-yNN pieces of 
that contribute to 7iV scattering at N 2 LO can be obtained from the 5th, 9th, 12th 
and 13th terms of Eq (3.8) and the last term of Eq. (3.9) in that paper. 

Finally, the relevant N 3 LO terms are given by 

C {3) = 2nN^ { [5/3 N + 5(3 v r 3 ] \g^ v - [(5a N + 5/3 N ) + (6a v + 6f3 v )r 3 ] v^} F» p F v p N + . . . , 

(15) 

where 5a^ (SPn) and 5a v (5f3 v ) are short-range contributions to the isoscalar and isovector 
electric (magnetic) polarizabilities of the nucleon, 5a p = (5o>n + 5a v ) /2 and 5a n = (^a^ — 
8a v )/2 (similarly for 5(3 P and 6{3 n ), which we seek to determine. These LECs are linear 
combinations of e 8 Q-e 9 4 of Fettes et al. [HE] . The fixed-coefficient pieces of e 8 9-e 94 are given 
in Table 6 of that paper. Other relevant terms in which have only fixed coefficients 
and are not written in Eq. (|THj) but can be found in Ref. fSJ are Om and Ons of Table 
5, as well as X 41 and X 53 of Table 8 and Y"n of Table 9. 
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2.2 Power counting 



Scattering amplitudes are computed from Feynman diagrams built from the effective 
Lagrangian. Observables are then power series in momenta and pion masses, with the 
non-analyticities mandated by perturbative unitarity. The generic form for an amplitude 
is shown in Eq. ([[J. 

Counting powers of Q is simple, except for one important subtlety The simple part is 
that each derivative or power of the pion mass at a vertex contributes Q to the order of 
a diagram, each loop integral contributes Q A , each delta function of four-momentum con- 
servation contributes Q~ A , and each pion propagator contributes Q~ 2 . The complication 
arises from the fact that in multi-nucleon systems there exist intermediate states that differ 
in energy from the initial and final states by only the kinetic energy of nucleons, which is 
0(Q 2 /M). Diagrams in which such intermediate states appear are called reducible, while 
ones containing only intermediate states that differ in energy from initial/final states by 
an amount of 0(Q) are called irreducible. 

First we consider only irreducible diagrams. By definition these are graphs in which 
the energies flowing through all internal lines are of 0(Q). A nucleon propagator then 
contributes Q~ x , and we find that a diagram with A nucleons, L loops, C separately 
connected pieces, and Vi vertices of type i contributes Q u , with [T3J[3T] 



Eq. (JJJ follows immediately from the assumption that the parameters in the effective 
Lagrangian scale as numbers of 0(1) times the appropriate power of A x . Hence the 
leading irreducible graphs are tree graphs (L = 0) with the maximum number of separately 
connected pieces (C = max), and with leading-order (A* = 0) vertices. Higher-order 
graphs (with larger powers of v) are perturbative corrections to this leading effect. In this 
way, a perturbative series is generated by increasing the number of loops L, decreasing 
the number of connected pieces C, and increasing the number of derivatives, pion masses 
and/or nucleon fields in interactions. 

Processes involving one or zero nucleons and any number of light particles with p ~ m n 
only receive contributions from irreducible diagrams, and there is good evidence that the 
resulting %PT power counting defined by Eq. (fT7)|) works in most kinematics (exceptions 
are possible, e.g., near threshold pion three-momenta may be smaller than the ~ m n 
assumed in deriving Eq. (j!6|l ) In processes involving more than one nucleon, the full 
amplitude consists of irreducible subdiagrams linked by reducible states. New divergences 
arise from these reducible intermediate states, which can potentially upset the scaling of 
the EFT parameters with A x . Weinberg ^3] has implicitly assumed that this is not the 
case, and thus that the power counting continues to hold for irreducible diagrams. 




(16) 
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It is now known that this assumption is not always correct [121 • m particular, consis- 
tency of renormalization in NN scattering seems to require that some pion-mass contri- 
butions should be demoted to higher orders than in Weinberg's original power counting. 
On the other hand, the error of not doing so is apparently numerically small [IS]- This 
explains the phenomenological successes of applications of Weinberg's proposal to few- 
nucleon systems |7| I41j. including to 7d scattering ^Hj and a number of other reactions 
involving external probes of deuterium . Here we continue to organize our calculation 
according to the power counting (|TB|) . but we will see some indications that the irreducible 
diagrams for Compton scattering on deuterium are not behaving exactly as Weinberg's 
straightforward chiral power counting for them would suggest. 

In an A-nucleon system which is not subject to any external probes, the sum of A- 
nucleon irreducible diagrams is, by definition, the potential V. To obtain the full two- 
nucleon Green's function G, these irreducible graphs are iterated using the free two- 
nucleon Green's function Gq, which can, of course, have a small energy denominator 
that does not obey the ;^PT power-counting: G = G + G VG. Making a spectral 
decomposition of G we can extract the piece corresponding to the deuteron pole. Since V 
is two-particle irreducible we may perform a chiral power counting on it as per Eq. 
Such a potential was constructed including terms up to v = 3 and with an explicit Delta 
field in Refs. jH]. Better fits have now been achieved using a potential without explicit 
Delta's at v — 3 |1q"1 14*6] . Recently the first calculation of V to v = 4 appeared, and a 
good fit to iVTV scattering data for lab energies below 290 MeV resulted jl?] . Here we use 
the wave function obtained from the NLO [y = 2) potential of Ref. [13], as that is the 
consistent choice given the chiral order at which we calculate 7d scattering. Since this is 
only an NLO potential the quality of the fit to NN scattering data is not nearly as good as 
that found in modern phenomenological potentials, or in Ref. j57j . In order to get a sense 
of how large the impact of higher-order terms in V might be, we have also computed 7d 
scattering using deuteron wave functions generated from a modern one-boson-exchange 
NN potential 0, and the Bonn OBEPQ @|. 

When we consider external probes with momenta of order Q, the sum of irreducible 
diagrams forms a kernel K for the process of interest, here the kernel i^ 77 for 'jNN — > 
jNN. Again, the power counting of Eq. (|T^|) applies to this object, since all two-nucleon 
propagators in it scale as Q~ l , as long as the photon energy uj ~ Q. Note also that in 
order to have an irreducible graph in which two photons are attached to different nucleons 
the two nucleons must be connected by some interaction, so that there is a mechanism 
to carry the energy of order Q from one nucleon to the other. Thus such graphs should 
be regarded as having one separately connected piece {C = 1), so that the dictates of 
four-momentum conservation can be respected. 

The full Green's function G 77 for the reaction is now found by multiplying the kernel 
by two-particle Green's functions in which the NN propagators with small energy 
denominators may appear, G 77 = GK J7 G. Extracting the piece of G 77 corresponding to 
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the deuteron pole at E = —B in the both the initial and final state we obtain the 7d 
T-matrix: 



Full consistency can be achieved by treating both the kernel for the interaction with 
external probe and the NN potential V in xPT and expanding both to the same order 
in v. Note that if is constructed in a gauge invariant way and is worked out to the 
same order in xPT as V then this calculation will be gauge invariant up to that order in 
%PT. It will not, however, be exactly gauge invariant. Contributions to K^ n of the form 
K^GK^ (with K 1 the two-nucleon-irreducible kernel for one photon interacting with the 
iViV system) only appear as part of i^ 77 order-by-order in the expansion in Q/A x . For 
gauge invariance to be exact this contribution from two-nucleon intermediate states must 
be included in full in G 77 . 

This is not a deep problem. In the very-low-energy region u ~ rn^/M, the power- 
counting of Eq. (fTT)|) does not apply to K ll7 since graphs in K 77 with two-nucleon inter- 
mediate states will have denominators of order Q 2 /M. There are thus two regimes for 
Compton scattering on nuclei, depending on the relation of the energy of the photonic 
probe, u, to the typical nuclear binding scale ~ m^/M. (Compare with the discussion 
of low-energy theorems for threshold pion photoproduction on nuclei in Ref . jH] , and of 
Compton scattering on the deuteron with perturbative pions in Ref. (SHI-) Here we are 
interested mostly in photon energies between 60 and 200 MeV, and so we regard ourselves 
as being firmly in the regime uj ~ m n 3> m^/M. In this regime it is valid to treat the 
contributions K^GK^ perturbatively in the Q expansion, i.e. we can use a chiral expan- 
sion for the two-nucleon, fully- interacting, G which appears in K^GK^. This results in 
a calculation that is completely gauge invariant up to the order in the chiral expansion 
to which we work. However, such a treatment fails once uj ~ m^./M, and so — as will 
become manifest below — our conclusions regarding 7d scattering at these lower energies 
should be viewed with some caution. In this very-low-energy domain it seems sensible to 
integrate out the pion, and then analyze mechanisms and count powers in the "pionless" 
EFT HH- 




(17) 




(18) 
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2.3 Notation 



A word of explanation about notation is required. We are interested here in Compton 
scattering on both the nucleon and the deuteron to next-to-next-to-leading order, since 
it is at that order that short- distance contributions to the polarizabilities first appear. 
These short-range jN interactions have index A = 3, see Eq. ([15)1 . On the other hand, 
the corresponding value of the index v defined in Eq. (116)1 depends on the target. 

For scattering on a single nucleon, A = 1 and C — 1 in Eq. Qlbjl. The first non- 
vanishing contribution to the T-matrix — the Thomson term — appears at v = 2, while 
the short- distance polarizability terms first show up at u = 4. One-nucleon scattering 
diagrams are traditionally labeled as Q u . The highest order we work to here is thus 
0(Q 4 ). 

For scattering on the deuteron, A = 2 and C — 1,2 in Eq. (|16|) . The contributions 
of the single-nucleon amplitude to K yi have C — 2. The single-nucleon Thomson term 
therefore has v — — 1, while the short-distance polarizability terms first appear at v — 
1. However, in order to make closer contact with previous work on 'jN scattering in 
%PT |1(J[ ITU IT2] , we follow the single-nucleon practice and refer to our v = 1 deuteron 
calculation as being 0(Q 4 ). Of course only relative orders are physically significant. 

A final warning. Our index A counts simultaneously the number of derivatives and 
nucleon fields. Sometimes the EFT Lagrangian is broken down into a purely-pionic La- 
grangian, a pion-nucleon Lagrangian, and so on, and the terms in each are further labeled 
not by A but by the number of derivatives. For example, our £(°) of Eq. (JJJ) is often 
written as C^l + C^ N + C^ N . These notational issues should be borne in mind when 
referring to Refs. [T^l ITS] . 

3 Compton scattering on the proton: theory 

In this section we consider Compton scattering on a single nucleon. This is required for 
the extraction of the proton polarizabilities from proton Compton scattering data, and 
also because the ;\FT 7 A" amplitude is an important ingredient in our calculation of 7d 
scattering. We work in the Coulomb gauge where the incoming and outgoing photon 
momenta k = (co,k) and k' = (to', k') and polarization vectors e and e' satisfy: 

k 2 = 0; k' 2 = 0; v ■ e = 0; v ■ e' = 0. (19) 

In either the Breit frame or the 7^ center-of-mass frame the Compton scattering 
amplitude can be written as: 

T lN = I Ate ' ■ e + A 2 e ' ■ke-k' + iA 3 a ■ (e ' x e) + iA 4 a ■ (k' x k) e' ■ e 
+iA 5 a ■ [(? ' xk)e-k'-(ex k') e ' ■ k] + iA 6 a ■ [(?' x k') e ■ k' - (e x k) e' ■ k] } , (20) 
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(d) (e) (f) 



Figure 1: Tree diagrams that contribute to Compton scattering in the e • v — gauge up 
to 0(Q A ). Small dots are vertices from larger dots are vertices from the sliced 
dots are vertices from and the sliced and diced dots are vertices from C^ 3 \ Crossed 
graphs and graphs which differ only in the ordering of vertices from £W an d are 
included in the calculation, but are not shown here. 



with Ai, i — 1 . . . 6 scalar functions of photon energy and scattering angle. 

In the relativistic theory the Compton scattering amplitude can be split into Born 
and non-Born — or "structure" — terms. The six polarizabilities a, (3, and 71-74 are the 
independent coefficients that arise when these structure terms are expanded in powers of 
uj up to uj 3 . For the Breit-frame amplitude such an cu-expansion gives (keeping terms up 
to 0(1/M 3 )): 



Ai(u;,0 

A 2 (lu,9 
A 3 (u,6 

A 4 (u,9 
A 5 (oj,9 
A 6 (lu,6 



Z 2 e 



2„2 



M +4^( (Z + K)2(1 + C ° S ' 



Z 2 )(l -cos6) uj' 



AM 3 
e 2 oj 



(a + (3 cos 6)u 2 + 0(lu a ) 
k(2Z + k)uj 2 cos 6 - 4irf3u 2 + 0(to A ) 



, )u2 (Z(Z + 2k)-(Z + k) 2 cosO) +Af 

+4770/(71 - (72 + 2 74 ) cos9) + 0(lu 5 ) 



e 2 uj 

~2M 2 

e 2 oj 



\Z + k) 2 + AI +4nu 3 ~? 2 + 0{uj' 



2M 2 
e 2 uj 



2M 2 



[Z + k) 2 + AI +47tcj 3 74 + 0(cj 5 ) 
Z{Z + k) + Af + 4tt^ 3 7 3 + 0{u h 



(21) 



with A^°-Aq° the contributions from the pion-pole graph, Fig. ^d). 

In heavy-baryon chiral perturbation theory the amplitudes A\-A G include the tree 
contributions shown in Fig. ^ and the pion-loop contributions of Fig. EJ At 0(Q 2 ) only 
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Figure 2: Diagrams which contribute to nucleon Compton scattering in the e ■ v = 
gauge at 3rd (i-iv) and 4th (a-s) order. Vertices are labeled as in Fig. ^ 

graph CJa) contributes, and this reproduces the spin- independent low-energy theorem for 
Compton scattering from a single nucleon — the Thomson limit: 

T,n = ~e- r'^f , (22) 

where Z is the nucleon charge in units of |e|. 

At 0(Q 3 ) the s-channel proton-pole diagram, Fig.^b), and its crossed u-channel part- 
ner, together with a 'jN seagull from (Fig.^c)) ensure that HB%PT recovers the Low, 
Gell-Mann and Goldberger low-energy theorem for spin-dependent Compton scattering 
|4Uj . At 0(Q 3 ) the t-channel pion-pole graph Fig. EJd) also enters; its contribution, which 
varies rapidly with energy, is often included in the definition of the spin polarizabilities 
and is the reason that the backward spin polarizability is so much larger in magnitude 
than the forward one. Pion-loop graphs enter first at 0(Q 3 ) — see graphs i-iv of Fig. |2J— 
and give energy- dependent contributions to the amplitude which include the well-known 
predictions for the spin-independent polarizabilities of Eq. (jSJ), as well as less- famous pre- 
dictions for 71-74 [TUj. The third-order loop amplitude of Refs. |§J H2] is obtained from 
diagrams llfi)-[2ti v ) — together with graphs related to them by crossing and alternative 
time-ordering of vertices — and is given in Appendix [X] As observed in Ref. [25j, the 
0(Q 3 ) xPT result for 7p scattering is in good agreement with the data at forward angles 
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up to quite high energies (~ 200 MeV), but fails to describe the backward-angle data 
already at around 100 MeV. 

At 0{Q i ) in HB^PT nucleon pole graphs ^e) and the fixed-coefficient piece of the 
seagull depicted in Fig. ^f) give further terms in the 1/M-expansion of the relativistic 
Born amplitude. There are also a number of fourth-order one-pion loop graphs (graphs a-r 
of Fig. I2J), which — with one exception — contribute only to the "structure" part of the 7A" 
amplitude. (The exception is graph |2k which also renormalizes k in the spin-dependent 
Born amplitude.) All of the graphs EM2t involve vertices from of Eq. (|T3~j) . In 
particular, the expressions for diagrams |2hH2t contain the LECs c±, c 2 , and c 3 , and in 
diagrams EM2k the nucleon anomalous magnetic moments enter. 

Divergences in graphs l2M2k, Hb) and Eh need to be renormalized by the coun- 

terterms appearing in Eq. (fTo^) . The divergent fourth-order loop contributions to the 
spin-independent polarizabilities — first derived in Ref. JT] — are as follows: 



loop • 4nf 2 
2 

e 



P] 



(l2M + M ~ 1} + (12/is + 2d + 4)T3 ) ~ l~° 2 ) K] 

* (24M (~ d2 + m + 23 + 24/isT3 ) " K 52 + {d ~ 4)(c3 ~ 2ci) )) A - 

(23) 



where d is the space-time dimension, 

A; = (d - 2)L(ji) + (1/Stt 2 ) log f^j , (24) 

H is the renormalization scale, and 

(d-4)L(fi) = (l/16vr 2 ) +0(d-4). (25) 
Now adding the counterterms from of Eq. ()15j) cancels the divergences, leaving: 

M* + PP = I ^(^((94 + 24(^ + l)r3)log^+79 + (16 + 12^)r 3 ) 



Sa N (fi) + 6p N (p) + r 3 {8a v (jj) + <%,(//)) ; 



4 ^> - v&n ( 4 (< 94 + 48 ^> lo « (t) + 51 + 24 "* rs 



2 fm \ 1 

-c 2 log ( -j\ - - (ca + 2c 3 - 4c x ) ) + 6p N (n) + TsSPM, 



(26) 



with g 2 = c 2 - 



16 



The scale dependence in the LECs 5a>N,v and 5(3n,v cancels the /^-dependence of the loop 
pieces. The 7 A?" seagulls proportional to these LECs (Fig. HIT)) represent contributions 
to a and (3 from mechanisms other than soft-pion loops, i.e. short- distance effects. As 
such they should, strictly speaking, be determined from data, since they contain physics 
that is outside the purview of the EFT. In contrast the spin polarizabilities, as well as all 
higher terms in the Taylor-series expansion of the Compton amplitude, are still predictive 
in xPT at this order — as long as the values of the Cj's are taken from some other process. 
The first short-distance effects in 71-74 occur at 0(Q 5 ). 

The full amplitudes at 0(Q 4 ) are too lengthy to be given here; they can be constructed 
from the information given in the appendix of Ref. [12] and the details given in Appendices 
|B]and0 We stress that in what follows it is always the full 0(Q 4 ) HB;\RT amplitude, 
and not a Taylor-series expansion such as (f2*T]) . which is employed. 

A well-known problem with the heavy-baryon approach is that at finite order the pion- 
production threshold comes at uj = m n . Furthermore beyond third order the amplitudes 
diverge at that point. The solution is to resum the terms which, taken together to all 
orders, would shift the threshold to the physical point uj^; in practice this involves writing 
the amplitudes as functions of uj/uitk instead of uj /m n , then Taylor-expanding to correct 
the error introduced to appropriate order. The resulting amplitudes are then finite and 
have the threshold in the correct place. This procedure was first introduced for forward 
scattering in Ref. ^0], and elaborated for non-forward scattering in Ref. ^2]- Since the 
proton data fitted spans the pion-production threshold we used this procedure for the 
extraction of the proton polarizabilities, but as the highest energy in the deuteron case is 
only 95 MeV, we did not do so in that case. 

Having obtained the amplitudes, the differential cross section in any frame is a kine- 
matic prefactor $ 2 times the invariant |T| 2 . Experimental data, though taken in the lab 
frame, are sometimes presented as c.m.-frame differential cross sections, so both prefactors 
are needed: 

•» = r w • *» = V-r- < 27 > 

Air 47T y's 

Finally, in terms of the Breit-frame Ai-A$, \T\ 2 is given by: 

|T| 2 = ±A 2 (l + cos 2 #) + ±A 2 (3-cos 2 #) 

+ sin 2 9 [AA 3 A 6 + (A 3 A 4 + 2A 3 A 5 - A X A 2 ) cosO] 

+ sin 2 6 \\A\ sin 2 9 + |A 2 (1 + cos 2 9) + A 2 5 {1 + 2 cos 2 9) + 3A 2 6 

+2A 6 (A 4 + 3A 5 ) cos9 + 2A±A 5 cos 2 9]. (28) 

(Above the nN threshold one should replace A 2 by \Ai\ 2 and AiAj by $l(A*Aj).) 
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4 Compton scattering on the proton: results 



The only unknowns in the 0{Q A ) ;\RT 7p amplitude are 5a p and 5j3 p , since the values 
of c\, and C3 can be determined from ttN [39 j (or NN |4L)j ) scattering data. If these 
counterterms are combined with the (scale-dependent) divergences they renormalize we 
can regard a p and f3 p (or, more precisely, the shifts of the polarizabilities from their 0(Q 3 ) 
values) as the two free parameters in the 0(Q A ) calculation. 

In Ref. ^2], these two free parameters were not fit directly. Rather, the input to that 
calculation included the Particle Data Group (PDG) values of the polarizabilities [HQ] . 
The resulting differential cross sections are in good agreement with the low-energy data. 

While these PDG values are derived from polarizabilities quoted in experimental pa- 
pers, the energies of the relevant experiments are, for the most part, too high for a direct 
fit to the polynomial form (j21j) . Instead the amplitudes these papers employ in order to 
extract a p and [3 P are based on dispersion relations, using plausible assumptions about 
their asymptotic form coupled with models for those amplitudes which do not converge. 
The difference a p — (3 P is a parameter in these fits; the sum is (usually) taken from the 
Baldin sum rule j^TJ. 

This approach implicitly builds in model-dependent assumptions about the behavior 
of the 7p amplitude beyond the low-energy regime. In what follows we adopt a different 
strategy, using x?T to determine the polarizabilities from only low-energy (u < 200 MeV) 
differential cross sections, so as to obtain truly model-independent results for a p and (5 p . 

In a recent paper the world data on proton Compton scattering was examined, and 
it was demonstrated that data from the 1950s and 1960s is compatible with the more 
modern data from 1974 on, and is useful in reducing errors |H2j- We have therefore used 
all the data in Table 2 of Ref. [52] , with the exception of the first three points for which 
the normalization error is not known. We have also included higher-energy points from 
the experiments of Refs. ^7j, and the TAPS data [T5] . 

All the experiments have an overall normalization error. We incorporate this by adding 
a piece to the usual \ 2 function: 



where dcr^ and Ajj are the value and statistical error of the ith observation from the 
jth experimental set, 8j is the fractional systematic error of set j, and Mj is an overall 
normalization for set j . The additional parameters Mj are to be optimized by minimizing 
the combined x 2 - This can be done analytically, leaving a \ 2 that is a function of a p and 
(5 P alone. Since we fit 13 7p data sets this reduces the number of degrees of freedom in 
the final minimization by 13. (For more details on this formalism see Ref. |52j.) 
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co max (MeV) 


Adata 


a p (10" 4 fm 3 ) 


fa (lO" 4 fm 3 ) 


X 7d.o.f. 


200 


144 


11.7 


3.40 


1.32 


180 


126 


12.1 


3.38 


1.18 


165 


109 


12.3 


3.37 


1.27 


140 


73 


12.0 


3.13 


1.26 



Table 1: Results of 0(Q 4 ) fits for a p and P P to world 7p data over various ranges u, \f\t\ < 
u max . In each case A^ata is the number of data points included in the fit. 

In Ref. ^2] it was found that while the xPT fit was very good at sufficiently low 
energies and forward angles, it failed for higher energies and backward angles, and did so 
in a way that could not be remedied by adjusting parameters. This suggests that at a 
relatively low momentum "new" physics is coming in — an observation that is consistent 
with the Delta- nucleon mass difference being only a factor of two larger than m^. Indeed, 
in both existing models and the effective field theory with explicit Delta-isobar degrees 
of freedom [SSI Ell Delta effects appear at backward angles already at uj ~ m n . (Their 
impact on the forward- angle data does not show up until higher energies are reached.) As 
a consequence in our extraction of a p and P p we impose a cut on both the photon energy 
and the momentum transfer, oj, yfi[ < uj max . 

We investigated the effect of varying cu max between 140 and 200 MeV — see Table [TJ 
As t^max was increased the central values of a p and (3 P changed very little; the 1-a curve 
tightened as more data were included, but the xVd.o.f of the best-fit point tended to 
rise. This did not seem to be due to a systematic worsening of the fit, but rather to the 
inclusion of rather noisy data. We have chosen to quote the central value corresponding 
to a cut-off of uj m &x — 180 MeV, for which x 2 /d.o.f was (locally) minimized at 133/113. 
This fit is shown in Fig. Altogether with the low-energy data. 

The best-fit point for the proton electric and magnetic polarizabilities is shown in Fig. 
IU together with the 1-a curve (Xmm + 2.28): 

a p = (12.1 ± l.l)loj x 10" 4 fm 3 , 

P P = (3.4 ± x 10" 4 fm 3 . (30) 

Statistical (1-a) errors are inside the brackets. An estimate of the theoretical error due 
to truncation of the expansion at 0(Q A ) is given outside the brackets. The result (|5Dj) 
is fully compatible with other extractions, although the central value of P P is somewhat 
phigher [Hilt IE3]. 

This theoretical error is assessed by varying the bound on which data are fit from 140 
MeV to 200 MeV. The size inferred thereby is consistent with estimates of the impact that 
0(Q 5 ) terms would have on the fit. In this context it is worth noting that although the 
0(Q 4 ) ;\RT predictions for the spin polarizabilities 71-74 are not in good agreement with 
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Figure 3: Results of the 0(Q A ) EFT best fit to the differential cross sections for Compton 
scattering on the proton at various angles, compared to the experimental data fTJ . 
The gray region is excluded from the fit. The magenta diamonds are Mainz data ^H]; the 
other symbols are explained in Ref. [T2*] . 



those obtained from other sources, in Ref. [12J it was shown that at the energies considered 
here the differential cross section is largely insensitive to their values. (Asymmetries in 
polarised Compton scattering would be far more sensitive.) 

The recent re-evaluation of the Baldin sum rule [HI] by Olmos de Leon et al. [IH] gives: 

a p + p p = (13.8 ± 0.4) x 10~ 4 fm 3 . (31) 

If one considers only statistical errors the results (J3Uj) and (J3TJ) are marginally inconsistent 
at the 1-cr level — see Fig. |3J Our values for a p and (3 P are consistent with the constraint 
()31j) once the potential impact of theoretical uncertainties is taken into account. 

Including the constraint (J31)) in our fit leads to values for a p and (3 P consistent with 
Eq. fl30(l . but with smaller statistical errors, namely: 

a p = (11.0 ± 0.5 ± 0.2)i°j x 10~ 4 fm 3 , 

p p = (2.8 ± 0.5 =F 0.2)±g;J x 10" 4 fm 3 . (32) 

In Eq. (J22I) we have left the systematic error unchanged, but have now included a second 
error inside the brackets, whose source is the error on the sum-rule evaluation (|31|). The 
smaller statistical errors are achieved at the expense of introducing certain assumptions 
about the high-energy behavior of the Compton amplitude into the result, via the use of 
the Baldin sum-rule result (pHj) . The result (}3*2j) is in excellent agreement with a recent 
determination of a p and [3 P which used an effective field theory with explicit Delta degrees 
of freedom to fit all 7p scattering data up to uj ~ 170 MeV p?5]. 
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Figure 4: Values (in 10 _4 fm 3 ) of the proton magnetic and electric polarizabilities ex- 
tracted from an 0(Q A ) EFT fit to low-energy cross-section data, together with the values 
allowed by the Baldin sum-rule constraint ()31|) (gray shaded area). 

5 Compton scattering on the deuteron: theory 

In this section we concentrate on real Compton scattering on an A = 2 system. 

The leading contributions to the irreducible 'yNN —>■ 'yNN kernel appear at 
v = — 1 — 0(Q 2 ) — from a tree-level diagram (L = 0) with two separately connected pieces 
(C = 2) and a vertex which is the two-photon seagull of A, = 1. This vertex arises 
from minimal coupling in the kinetic nucleon terms in of Eq. (fT3*|) — see Fig. Efa). 
It therefore represents Compton scattering on the NN system where the neutron is a 
spectator while the photon scatters from the proton via the Thomson amplitude (|22j) . 

At orders v > there are two classes of corrections: corrections to the one-nucleon 
amplitude where the second nucleon is just a spectator, and corrections that involve both 
nucleons. Shown in Fig. I^b-e) are the one-nucleon corrections at v — 0, while Fig. |U] 
displays two-nucleon corrections of the same order. Note that the one-nucleon-loop and 
two-nucleon diagrams involve only Aj = interactions in diagrams with L — 1 or C — 1. 

At this, or indeed any, order the irreducible ^NN — > ^NN kernel can be separated 
into one- and two-body pieces. Doing this in the jNN cm. frame, with the kinematics 
shown in Fig. [7| we obtain: 

Kff^tf, k'-p, k) = T] d N c - m -(k'; k) 5® ( p '-p- lq) + T 2 J N (p>, k';p, k); (33) 

with q = k — k' the momentum transfer of the Compton scattering process. The one- 
nucleon piece T^ c ' m '(k'; k) is related to the amplitude discussed in Sec. El by a boost 
from the Breit frame to the 7d cm. frame. Details of the boosting procedure are given 
in the Appendices. 
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Figure 5: Characteristic one-body interactions that contribute to Compton scattering 
on the deuteron at order Q 2 (a) and at order Q 3 (b-e) (in the Coulomb gauge). Crossed 
graphs are not shown and not all loop topologies are shown. The sliced blobs are vertex 
insertions from C^ 1 '. The sliced and diced blobs are vertex insertions from C^ 2 \ 




Figure 6: Two-body interactions which contribute to Compton scattering on the deuteron 
at order Q 3 [y = 0) in Coulomb gauge. Graphs which differ only by interchange of the 
nucleons are not shown. 
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Figure 7: Kinematics for the 70! process in the 70I center-of-mass frame. The sliced circle 
represents K^, the irreducible ^NN — > ^yNN kernel. 



Meanwhile T$ N was computed to third order in small quantities in Ref. ^Hl- At fourth 
order the difference of initial and final nucleon kinetic energies should be included in the 
pion propagators appearing in Fig. |U1 The result is then: 



with: 



tw = ^ +(1 „ 2) 

2[c<j 2 — m 2 — (p — p' + k + ) 2 ] 
t (b) = e-e' a 1 -{p-p' -\q)a 2 -(p-p' + \q) | ^ 2 
2[(p — p' - |<f) 2 + m 2 ][{p — p' + \q) 2 + m 2 } 

t ( c ) = g ' ■ (g - g + g ■ g g ■ (P ~ g + |g) | ^ ^ ^ 

[u; 2 - m 2 - (p* — p' + A;+) 2 ][(p - p' + |g) 2 + m 2 ] 

^) = 1 (P ~ g + jgO ^ • (P - g ! - ¥ ))^ • + (i <_> 2) 



£ (e) 



p 2 - m 2 - (p — p' + £; + ) 2 ][(p - p' - \q) 2 + m 2 ] 
a 1 ■ (p- p' - \q) a 2 ■ (p '— p' + ~<f ) 
— p'— \q) 2 + m%\[(p - P 1 + |<?) 2 + " m2 r ] 



x ^jgz^±Ml£^g^±M) + (1 „ 2) , (35 ) 
[uj 2 — m 2 — (p — p' + /c+) 2 ] 

where fc+ = + fe') is the average of incoming and outgoing photon three-momenta, the 
notation (1 «-> 2) indicates the inclusion of a second term in which p — > — p, p' — > — p', 
a 1 <-» (j 2 , and: 

— # — # 

ft 2 = „ 2 + ^'' fc/ ~ jr ' fc . (36) 

In computing the "retardation" effects which modify modify u — > uj we use the fact 
that — because we are considering an elastic process — the two-nucleon initial and final 
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states are equally distant from being on-shell. Thus the results (J33J), (pjj) were derived 
using \p | 2 = \p'\ 2 , which simplifies the expression for u. Note that for virtual momenta p 
and p ' of order M the retardation effects included here will have unphysical consequences, 
and thus Eq. (jHHj) is only valid if \p\, \p'\ -C M. 

Upon including these retardation effects in our calculation we found them to be nu- 
merically very small. Hence in the computations of the 7d differential cross section which 
we report on below they were not considered. 

Since neither the expression for X^vato(<2 3 ) nor ^ ne 7^ 0{Q Z ) amplitude contain any 
free parameters, xPT makes predictions for 7d scattering at this order. These predictions 
were generated in Ref. [T3] and compared to the data of Ref. [TP] , with encouraging results, 
especially at = 69 MeV. The 0(Q 3 ) predictions (made on the basis of Ref. |T3]) for 
the Lund experiments at uj pd 55 MeV and u ~ 66 MeV turned out to be in excellent 
agreement with that data too. (See curves below.) 

At 0(Q A ) there are many one-nucleon diagrams in the ^NN kernel, which are easily 
obtained from those in Figs. ^ and 121 These were computed in the Breit frame in Ref. [T2]. 
The extra pieces of the 7 A" amplitude which are obtained in the 7d center-of-mass frame 
as compared to the Breit-frame result are discussed in Appendices lAlfDl 

Two-nucleon diagrams at 0(Q A ) in Coulomb gauge are depicted in Figs. El and El These 
two-nucleon v — 1 diagrams have C = 1 with one insertion from C^ l \ 

Although they are nominally of 0(Q 4 ), all of the graphs in Fig. [HJgive zero contribution 
to the elastic 7d amplitude at this order. Diagrams (a), (a2), (d), and (d2) do not affect 
this process because i^ 77 is of isovector character. Meanwhile graphs (b), (b2), (c), (e), 
and (e2) all involve the tiNN vertex from C^ 1 '. With the choice v = (1,0) the Feynman 
rule for this vertex reads: 

^J-J-iP + P')^, (37) 

where p (p') is the nucleon momentum before (after) the emission of the pion, and q is the 
emitted pion's four-momentum. In the context of the jNN kernel being discussed here 
the emitted pion is a "potential" pion jSI], and as such has go ~ p 2 /M. Consequently 
in the kinematics we are considering all five of these graphs are suppressed to 0(Q 5 ). 
Lastly, graphs (k) and (1) represent time-orderings where the photon interacts with a 
single nucleon while the pion is "in-flight". If, as is the case here, an instantaneous 
potential is employed to generate the deuteron wave function then consistency between 
the potential V and the 'jNN — > ^yNN kernel requires that graphs (k) and (1) not be 
included in K 77 . 

The non-zero two-body graphs at 0(Q 4 ) are depicted in Fig. El The (7d cm. frame) 
expressions for these pieces of the irreducible ^yNN — > ^yNN kernel sum to: 
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Figure 8: Some of the two-body graphs that are nominally 0(Q 4 ) for 70I scattering in 
Coulomb gauge. Ultimately all of the graphs shown here are either identically zero or 
suppressed to 0(Q 5 ) by kinematic effects, as explained in the text. Non-zero two-body 
contributions to 7d scattering at 0(Q 4 ) are shown in Fig. El 
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Figure 9: The two-body diagrams contributing to 7d scattering in Coulomb gauge at 
0(Q 4 ). 
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with: 



tif+9) = - 5fF ^tUt 21 {€^-(p + p / -fc + ) 

2[(p-p + ^qY + mi\ 
x -(P-P'+ \q)o 2 -e + a 2 -(p-p' + )a x • e] 

+ or- (e' x fcOO? 1 + <? 2 ) • (p - p" + if)} 

tih+92) = ou- ^tWt 21 {e-(p + p'-fc+) 
2 i(P ~P + f<?) + m 7rJ 

x [a 1 -(p-p' + \q)a 2 ■ e' + a 2 -(p-p' + |f ^ ■ e'\ 

-i(l + /^)e" • (ex ^(a 1 + a 2 ) • (p - p' + if)} 

^h-H+ha+^a) = _ 1 ? / . (p> + ^/ + £ + ) 

2[u; 2 — m 2 — (p — p' + A:+) 2 ] 

x • (p - p' + £+)^ 2 ■ e + a 2 ■ (p - p' + fc + )o !1 ■ e] 
-z(l + n v )(p- p' + • (e' x £') (a 1 + a 2 ) • e 
— e • (p + p' — k + )[a L ■ (p — p' + fc+)<? 2 ■ e' 

+a 2 -(p-p' + k + )a 1 -e'] 
+i(l + K„)(p - p" + • (ex ^(a 1 + a 2 ) • e'} 
t (j +m ) _ 1 1 

[u 2 -m 2 n -{p-p' + k + y\ K + (p - p' - 1?) 2 ] 

xe- (p-p' + £;+){e" • (p + p" + 

x [a 1 • (p-p' + k + )a 2 ■ (p - p" - i<f) 

+a 2 ■ (p-p' + k + )a x ■ (p-p' - |f)] 
+ K v )(e' x £') • (p - p" + fc + )(a 1 + a 2 ) • (p - p" - iff)} 

^(n+o) _ 



[cu 2 - ml - (p- p' + k + f] [ml + (p - p' + |f) 2 ] 
xe' • ip — p' + k + ){e ■ (p + p' — k+) 
x ' (p-p" + £+)^ 2 • (p - p" + \q) 

+a 2 ■ (p-p' + k + )a 1 ■ (p-p' + \q)} 
-i(l + n v )(p- p' + k+)(e x ^(a 1 + a 2 ) • (p - p" + if)}. (39) 



(In this case the contribution from diagrams related to those in Fig. |§] by interchange of 
the two nucleons has been explicitly included in the expressions.) 

The vertices from the next-to-leading-order nN Lagrangian which appear in these non- 
zero two-body diagrams all involve El and Ml couplings of the photon to the nucleons. 
Thus the only parameters appearing are the proton charge and the proton and neutron 
anomalous magnetic moments. There are no free parameters in the two-body 0(Q 4 ) con- 
tribution. This is good, as these mechanisms are of the same order as the counterterms 
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5aN and 5(3n which we are trying to fit. Extracting the polarizabilities from deuteron 
Compton data is much more straightforward if the 0{Q i ) two-body currents can be ex- 
pressed in terms of known parameters. 

To calculate the amplitude for 7d scattering, T jd , we first sandwich between 
deuteron wave functions and use the decomposition of Eq. (jHSJ). This yields (in the 
7d center-of-mass frame): 

T M'x'Mx(k', k) = J ^* M , (p + \q) T^;™ik', k) ^ M (p) 



+ 



r M >(p') T^ Nyx (k', k) ^ M (p) (40) 



where M (M') is the initial (final) deuteron spin state, and A (A') is the initial (final) 
photon polarization state, and k (k') the initial (final) photon three- momentum, which 
are constrained to \k\ = \k'\ = uo. 

The laboratory differential cross section can then be evaluated directly: 



da 



dQ L 16 



1 ( E i\ 1 \^ \ T ~,d 2 
j7r 2 U 6 2^ ^M'X'Mxl 
v 77 M'X'MX 



MA I ' 



(41) 



where E 1 is the initial photon energy in the laboratory frame, and is related to lu, the 
photon energy in the 7d center-of-mass frame, via: 

u = - El ; (42) 
y/1 + 2E 1 /M d 

and E' is the final photon energy in the laboratory frame: 



To evaluate the center-of-mass-frame differential cross section we employ: 

da _ M d 
dnZ = 167r 2 (M d + 2E 7 ) 6 



M'X'MX 



6 Compton scattering on the deuteron: results 

In this section we present our results for the cross section for Compton scattering on 
the deuteron including the one-nucleon and two-nucleon mechanisms described above. 
Our calculation therefore represents the full v — 1, or 0(Q 4 ), result for Compton 
scattering on the deuteron. We present these 0(Q 4 ) results, examine how far they depend 
on the choice of deuteron wave function, and discuss our calculation's breakdown as the 



28 




a 

-a 





20 
15 
10 

5 
0. 




90 



180 




O E=69 MeV 
O E=66 MeV 



90 

6 (deg.) 

cm v ° ' 



180 





i 

20 
15 
10 
5 





90 180 

9 m (deg.) 




6 (deg.) 

cm v ° 7 



Figure 10: Results of the 0(Q 2 ), 0(Q 3 ) and 0(Q 4 ) calculations of (lab and cm. as ap- 
propriate) differential cross sections for Compton scattering on deuterium at four different 
lab photon energies: 49, 55, 67, and 94.5 MeV. The data are from Illinois ^H] (circles), 
Lund [2T] (diamonds) and SAL (squares) [20] ■ The solid line is the 0(Q 4 ) calculation with 



a N = 13.6 x 10" 4 fm , fl N = 0.1 x 10~ 4 fin 3 (Fit I of Table EJ). The dot-dashed line is the 
(parameter-free) 0(Q 3 ) calculation, and the dotted line is the result of the leading-order 
[0(Q 2 )] calculation. 



photon energy is reduced to energies comparable to the nuclear binding scale, m^/M. We 
then present results for different fits to the database of recent 7d experiments. Because 
the deuteron is an isoscalar, the only free parameters are the isoscalar combinations 
5aN = (Sotp + 5a n )/2 and 8(3n = (5{3 P + 5f3 n )/2. In fact, of all the additional terms which 
appear in the 7d calculation when we go from 0(Q 3 ) to 0(Q A ), only these LECs affect 
the cross section significantly. We conclude the section by comparing our results with 
numbers for neutron polarizabilities obtained via other techniques. 

6.1 Results at 0{Q A ) and convergence of xPT 

Our results for Compton scattering on deuterium at 0(Q 4 ) are presented in Fig.HU! where 
the solid line represents the calculations at photon energies of 49 MeV (lab), and 55, 66, 
and 95 MeV (cm.). Also shown are the results of lower-order calculations at the same 
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energies, as well as data from Refs. ^11201 [H]- The wave function employed was the NLO 
%PT wave function of Ref . > with the cutoff A chosen to be 600 MeV. As emphasized 
above, the only free parameters in this calculation are the isoscalar polarizabilities oln 
and Pn- The 0{Q A ) plots in the figure were generated with the values for and (3^ 
obtained in Fit I of Table El below. However, before discussing this, and our other fits, in 
detail, we want to use these results to address a number of theoretical issues associated 
with our calculation. 

One of the strengths of effective field theory is that it establishes a perturbative ex- 
pansion for S-matrix elements — see Eq. (0). Having employed an EFT to compute 7d 
scattering it behooves us to ask whether the perturbation expansion is behaving as ex- 
pected. To this end in Fig. [TUlwe also plot the 7d cross section at leading (0(Q 2 ), LO, 
dotted line), and next-to-leading (0(Q 3 ), NLO, dot-dashed line) order. The correction 
from LO to NLO is sizable, but is consistent with an expansion parameter M £_ M ■ The 
effect in going from NLO to N 2 LO is surprisingly small, perhaps because two-body effects 
that enter are controlled by, at worst, u/M, which is significantly smaller than the nom- 
inal expansion parameter to /{Ma — M). Generically, the largest effect at N 2 LO comes 
from the shifts of the polarizabilities from their 0(Q 3 ) values. Since this effect goes as uo 2 
the shift in the 0(Q 4 ) result compared to that at 0(Q 3 ) is much more noticeable at 95 
MeV than at any of the lower energies. 

All of this is encouraging for our attempts to extract and (3n from the coherent 
Compton cross section. On the other hand, a skeptical view of the results in Fig. El 
would be that another order must be calculated in order to ensure convergence. We will 
return to this point below, but we note that the as-yet-uncomputed fifth-order x?T result 
for the 7 A" amplitude is a key element in any such calculation. The 0(Q 5 ) result for 7d 
scattering therefore requires a two-loop computation in the single-nucleon sector. 



6.2 Wave-function dependence 

We have computed the 7AW — > 7AW kernel to two orders beyond leading order. For 
computation of the matrix element that enters observables the NLO wave functions of 
Ref. jlH] are therefore a consistent choice since they too include effects of 0{Q 2 ) beyond 
leading. (They are "next-to- leading order" for AW scattering, since the 0(Q) corrections 
to AW vanish in the parity-conserving part of the AW potential.) 

Ref. [13] employed a cutoff in its NLO calculation: a cutoff that was varied between 
500 and 600 MeV. The results of our 7d computation do not depend markedly on this 
cutoff. However, our results are rather sensitive to which AW potential is chosen in 
order to generate the deuteron wave function. In particular, the Nijm93 AW potential- 
model [5| gives a deuteron wave function which, when used in Eq. (|40p . results in cross 
sections somewhat larger than those found with the wave function of Ref. (3^]. These 
two results are extremal, in the sense that other wave functions (e.g. the Bonn OBEPQ 
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Figure 11: 0(Q 4 ) calculation of cm. differential cross section for 70! scattering with 
uj = 66 MeV. Data from Illinois and Lund, as in Fig.^3 The dashed line is the (parameter- 
free) 0(Q 3 ) calculation with the NLO wave function. In this curve the polarizabilities 
have their 0(Q 3 ) values. All other curves were obtained using the 0(Q A ) ^NN — > 'jNN 
kernel, with = 13.6 x 10~ 4 fm 3 and (3^ = 0.1 x 10~ 4 fm 3 . The solid red line is the 
result with the NLO xPT wave function of Ref. [4*H] . and the short-dashed green, dot- 
dashed blue and dotted orange lines use, respectively, the Nijm93, the Bonn OBEPQ, and 
R = 1.5 fm + OPEP [31] wave functions. In all cases the cutoff A is chosen to be 600 
MeV. 



wave function and the simple one-pion-exchange plus square-well wave functions of 
Ref. [S3]) generate cross sections which fall in between them. 

In Fig. ^2 we present a selection of these results for the case uj = 66 MeV. Different 
wave functions yield cross section predictions which vary by about 10-20%, with greater 
variation at backward angles. The variation with choice of wave function decreases slightly 
with energy. This wave-function dependence comes almost entirely from the matrix ele- 
ments of the two-body operators. The impulse-approximation piece of the cross section 
probes deuteron structure at momentum transfers of 200 MeV or less, and all of these 
wave functions give very similar results in that domain. 

Our syNN — > ^yNN kernel is valid only at low momenta. Thus, in all results presented 
here we have introduced a cutoff A on the momenta p and p f in the six-dimensional integral 
of Eq. (|4T)|) . The cutoff function employed is: 

f(k)=exp(-k i /A i ), (45) 
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Figure 12: One Feynman diagram which is of a higher order in the chiral expansion than 
we consider here, but must be included in the calculation of the 7d amplitude if T ld is to 
be exactly gauge invariant. The shaded blobs represent deuteron vertex functions. The 
small dots are vertices from the leading-order ttN Lagrangian, C^°\ while the larger ones 
are vertices from 

in imitation of the wave-function calculation of Ref . [13] . For all our results we chose A = 
600 MeV, so as to be consistent with the A = 600 MeV wave function of Ref. Provided 
A > 600 MeV, inserting the cutoff function in the matrix element evaluation changes the 
results with the NLO ;^PT wave function by less than 2%. It does, however, eliminate some 
high- momentum strength in the matrix elements evaluated with the Nijm93 wave function, 
thereby reducing the cross section. Not including the cutoff function (J45)) increases the 
spread of cross sections in Fig.^^to 30% or more at 66 MeV. Furthermore, without such a 
cutoff, wave functions from potentials with deeply-bound states, such as the NNLO wave 
function of Ref. |45j . generate cross sections 10-100 times larger than those displayed 
here. With the cutoff in place, and A = 600 MeV, the NNLO wave function of Ref. jlH] 
gives a result that falls within the band defined by the NLO wave function and the 
Nijm93 wave function. 

2 

6.3 Modifying the power counting for to ~ ^ 

No matter which wave function is chosen our 0(Q 3 ) and 0(Q A ) calculations overestimate 
the 49 MeV cross section, as measured at Illinois by Lucas [T^]. This would seem to be 
connected to the fact that at photon energies to ~ the power counting employed here 
breaks down and the contribution from two-nucleon-intermediate states must be included 
in full in the calculation. As things stand diagrams such as Fig. ^Jonly occur in the chiral 
expansion of the jNN kernel at 0(Q 5 ) and beyond, and without the inclusion of graphs 
like this one T 7 ^ is not gauge invariant. Thus the power counting we have used so far for 
7d does not recover the Thomson-limit amplitude for Compton scattering on deuterium 
as uj — > 0. 

These difficulties stand in contrast to the EFT(/) calculations of Refs. j2H], which re- 
cover the 7d Thomson limit and also reproduce the 49 MeV Illionis data. In the nomen- 
clature of Refs. EH] our calculation of 7d scattering is valid in "Regime II": the 
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Figure 13: Leading diagrams in Regime I for 7d scattering in EFT(jf). 



kinematic domain where u is of order the deuteron binding momentum: uj ~ 7 = V 'MB. 
The Thomson limit for 7d scattering emerges in a different regime: "Regime I" , which 
corresponds to photon energies of order the deuteron binding energy, i.e. uj ~ B. 

We now repeat the argument of Ref. [3U] m order to shed some light on how the 7d 
Thomson limit emerges in Regime I. There the leading diagrams for Compton scattering 
on deuterium are shown in Fig. EH The expressions for the operators to be sandwiched 
between deuteron wave functions can be written (in the 7d cm. frame): 



6 {a) 

0<9 



M 6 e 



M 2 ' ' P UJ + J^-B- p - e ' P 

w ~ 2M d M 



-,2 



1 



{c) = -tt^c-P re'-PI (46) 

M 2 F ^ p . (P-2fc+) 2 p 2 F ' v ; 

2M d 2M 2M 



where we have omitted the Ml piece of the jNN vertex, as well as the portions of the 
El vertex proportional to k, the photon momentum, since they are irrelevant to what 
follows. Note that nucleon recoil is included in the intermediate-state propagators here, 
i.e. they are not the standard heavy-baryon propagators, which would give only 1/uj and 
— 1/uj respectively. 

Evaluation of the operators O^-O^ between "zero-range" deuteron wave functions 

m = \/f (47) 

yields the result [30J: 

e 2 [4 7 f\q\\ 2 7 4 f / Mco . \ 3/2 / M^\ 3/2 11^ _ 

(48) 
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with \q \ = \f2u>(l — cos(0)). (Pieces of the graphs in Fig. IT31 down by a relative factor of 
uj/M must be omitted in order to obtain (|48|l .) Taking the limit u —>■ here yields the 
Thomson limit for Compton scattering on the deuteron: 



e 2 



T 7d (u = 0) = -— e" • e, (49) 

up to corrections due to deuteron binding, which are very small, being of order B/(2M). 

An important feature of this argument is that in EFT(^) an exact accounting of effects 
due to the scale \[Mxo can be performed. They produce the second term in Eq. pgj). As 
discussed in Ref. jTH], such contributions are higher order in our power counting, since 
if lu ~ m n they involve high relative momenta inside deuterium, which means that the 
deuteron wave function only enters the integral in regions where it is very small. Thus 
our calculation includes only the first term in Eq. (J48|) . However, as u — > effects from 
the scale a/ Muj become important, since the deuteron wave function is no longer being 
evaluated at a high scale. So, a power counting which regards only diagram (a) as the 
leading-order result, as ours does, is not correct in this very- low-energy regime. Indeed, 
with diagram (a) alone included, the amplitude (f4*HJ) is not recovered: a threshold cross 
section that is a factor of M^/M too large results. This is the source of our difficulty in 
reproducing lower-energy deuteron Compton data. 

In the absence of a full accounting of the effects of this scale in we can only 

2 

treat more carefully those diagrams which seem to become enhanced for uj ~ ^5-. Of 
these diagrams the argument of Ref. [SO] indicates that the most important are lT3T b) and 
113( c). Provided we are in the regime uj ^> B, we can evaluate these diagrams keeping 
only the contributions from the pole at \p \ 2 = Muj. Neglecting the effects of analytic 
structure at p ~ 7 and p ~ uj , which are suppressed by B/u and u/M respectively, we 
get, forintb): 

T s - Channe i po i e = 2 J^j M J J^k ^* (p ~ 5^') ?'-P S (\P\- ^Mu) t-P i> (p~ l S ; 

(50) 

If we assume, for the time being, that has only S-wave components, and again use the 
hierarchy y/uM ^> uj = \k\ = \k'\ then: 
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T s - Channe i po i e = (Muf 2 |v> (>/Mu) J t'-e. (51) 

Note that in the strict HBxPT expansion the imaginary part of the 7d amplitude is 
zero to all orders, provided that uj < m n . However, Eq. (joTj) provides a leading effect in 
Im(T) once we resum the recoil effects in the baryon propagator. Of course, a number 
of other mechanisms contribute to Im(T), but Eq. (J51|) does allow us to estimate its size. 
It indicates that the imaginary part has little effect on elastic 7d cross sections at the 
energies considered in this work. 
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An expression similar to Eq. ([51)1 . but with a result that is ultimately purely real, exists 
for diagram IT37c): 

2 e' ■ e. (52) 
Inserting the Fourier transform of Eq. (|47|) into Eqs. (Jo"Tj) and ()52|) we recover a result 

2 

that agrees with Eq. P%J) . as long as w ^> B. More generally, if a; ~ and so ^ is 
not small, we see that IT37 b) and (c) are numerically as important as the Thomson term 
da) 2 . 

And yet including diagrams IT3T b) and (c) at the same order as the Thomson term 
ITST a) violates the strict HB^PT power counting. It amounts to resuming an infinite set 
of recoil corrections for the two-nucleon propagator, all of which are higher-order in the 
original power counting. Below we refer to this as the "very-low-energy resummation" . 
To assess the impact of this resummation which enhances diagrams I13f b) and (c) in the 
region u ~ m\/M we evaluated these graphs numerically using the full expression (J5U|) . 
together with its crossed counterpart. Diagram ITBT c) proves to be the more important of 
the two, since it can interfere destructively with the result for diagram I13f a). provided 
that \il>(V Muj)\ 2 does not completely suppress its contribution. Ultimately it reduces the 
cross sections substantially at 49 MeV, but has little impact on them at 95 (and even 69) 
MeV [T3] . This last result is not surprising, in that at these higher energies T s _ channe i voie 
is suppressed by factors of the deuteron wave function at momenta of order yjMm n . 

It is clearly necessary to do this very-low-energy resummation when the deuteron is 
probed at very low photon energies: the baryons can no longer be treated as static for 

2 

lu < -rp, the nuclear binding scale, and the Thomson limit (}4T?|) will only be recovered if 
recoil terms for the baryons are included in the calculation. However, arbitrarily includ- 
ing just [UR b) and (c) ignores the possibility that in this very-low-energy regime other 
diagrams may also be enhanced compared to their u ~ m n importance. 

In fact, as emphasized in Ref. [IB], getting low-energy theorems correct can be quite 
complicated in a theory where the photon scatters from the deuteron's constituent nucle- 
ons and pions, just as recovering the Thomson limit for the proton can be quite difficult 
in a constituent quark model of proton structure [2j. As pointed out at the end of Section 
12.21 for the EFT discussed here in the very-low energy region the breakdown into reducible 
and irreducible diagrams changes: is no longer irreducible, and one needs to account 
for the appearance of the full two-nucleon propagator G between photon absorption and 
emission. Working this out in detail remains an important challenge. Although it bears 

2 

pointing out that once u < EFT(^) may be a more efficient way to calculate Compton 
scattering on deuterium, since at these low energies it is not necessary to include pions 
explicitly in the theory in order to achieve a good description of the data |2"9*] . 

2 The role of similar contributions from on-shell intermediate-state nucleons in neutral-pion photopro- 
duction on the deuteron has been emphasized by Wilhelm |55| . 
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6.4 Fits to deuteron Compton data 

Over the past decade three experimental groups have measured elastic 7d scattering. The 
data set includes: 

1 . The pioneering measurements by Lucas at Illinois, who took four data points at 49 
MeV and two at 69 MeV ^H]- The four 49 MeV points were taken in two different 
runs, and will, in principle, have uncorrelated systematic errors, therefore we choose 
to class these four points as from two different experiments. 

2. Five measurements taken at SAL at photon energies ranging from 84-105 MeV [2*U] . 
Data over this energy range was consolidated into one energy bin, and we quote it 
at a lab photon energy of 94.5 MeV. 

3. The very recent data set from MAX-Lab from Lund, which includes 9 points at 
approximately 55 MeV, and another 9 measurements at approximately 66 MeV [21] • 

In Table E] we list all of the data, including systematic and statistical errors. Note that 
in the table the angles and cross sections are in the center-of-mass frame for the SAL 
and MAX-Lab data, and in the lab frame for the Lucas data. Note also that since the 
higher-energy Lund data was taken at almost the same energy as the 69 MeV Illionis 
data in what follows we show them on the same plot, and compare them to calculations 
at Ery = 67 MeV. In the fits by which and (3^ were determined the energies listed in 
Table 121 were used to calculate the values for {da / dVL) t ^ eory . 

These fits seek to minimize the x 2 defined in Eq. (|29|h where systematic errors are 
accounted for via a floating normalization. (In the case of the world 7d data the index 
j of Eq. (}2*9"j) runs from 1 to 5.) Minimizing this \ 2 employing the NLO xPT deuteron 
wave function of Ref. jlHl we find a x 2 P er degree of freedom of 2.36, with a best-fit result 
of a N = 13.6 x 1(T 4 fm 3 and (3 N = 0.1 x 10^ 4 fm 3 . This fit, already shown in Fig. HOI 
has a rather high \ 2 ■ Examining the cross sections resulting from this calculation shows 
that the main contributions to the \ 2 come from the two backward-angle SAL points, and 
from the fact that the calculation severely over-predicts the 49 MeV Illinois data. Both 
of these problems are associated with the regime of validity of our calculation, which is 
valid for u ~ m n , but breaks down both at lower energies ~ rn^/M, where some of the 
Illinois data were taken, and for values of \f\t\ of order the Delta-nucleon mass difference, 
where the last two SAL points occur. 

In order to deal with this difficulty at higher \/]t\ we followed the same procedure as 
employed above for the proton data and eliminated points with lu, y\t\ > 160 MeV from 
the minimization. This excludes the last two SAL points from the fit, because our theory 
apparently does not include all the physics necessary to describe these points accurately. 
Refitting then reduces the x 2 /d.o.f. to 1.95, with best-fit values of ct^ and [3n equal to 
15.4 x 10~ 4 fm 3 and -2.3 x 10~ 4 fm 3 respectively. 
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54.9 


136.3 


18.4 


1.7 


1.6 




54.6 


136.8 


17.2 


2.0 


1.4 




55.9 


137.3 


21.0 


3.2 


2.2 




65.6 


44.4 


18.6 


2.4 


1.4 




65.3 


44.9 


16.0 


2.8 


1.4 




67.0 


50.3 


15.2 


1.8 


1.2 




65.3 


125.3 


15.3 


1.3 


1.4 




65.6 


127.8 


16.0 


1.2 


1.1 




67.0 


131.8 


14.2 


1.5 


1.0 




65.6 


136.5 


15.1 


1.7 


1.3 




65.3 


136.8 


12.6 


1.7 


1.8 




67.0 


137.5 


15.0 


2.7 


1.2 



Table 2: Deuteron Compton data sets from Illinois [THj . SAL j20], and Lund |2*T| . Cross 
sections and angles are quoted in the lab frame for the Illinois data and the cm. frame 
for the SAL and Lund data. Statistical and systematic errors are reported in nb/sr. 
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Figure 14: Results of the 0{Q i ) best fit (including very- low-energy resummation) to 
the (lab and cm. as appropriate) differential cross sections for Compton scattering on 
deuterium at four different lab photon energies: 49, 55, 67, and 94.5 MeV. The data are 
from Illinois JH] (circles), Lund [21] (diamonds) and SAL (squares) [20 . The error bars 
represent the quoted statistical (only) uncertainties of these measurements. The solid line 
is the 0(Q 4 ) calculation with a N = 13.0 x 10~ 4 fm 3 , p N = -1.8 x 10~ 4 fm 3 (Fit III). The 
gray area is the region excluded from the fit (a>, y/\t\ > 160 MeV). The dot-dashed line is 
the (parameter-free) 0(Q 3 ) calculation. 



While the restriction on the kinematic range of the fit affects the central values of 
and j3]\f it does not reduce the \ 2 markedly, because the agreement with the 49 MeV data 
is still poor. In contrast, implementing the very-low-energy resummation (Sec. 16. 3J) and 
then refitting produces a marked improvement in our description of the data — especially 
at energies around 50 MeV. The x 2 /d.o.f is reduced to 1.33, and the best-fit values for 
nucleon electric and magnetic polarizabilities are now 13.0 and —1.8 respectively. This fit 
(Fit III) is presented in Fig. and is clearly quite good; especially when one considers 
that the overall normalization of each set can be adjusted within their quoted systematic 
uncertainty. Omitting the 49 MeV Illinois data from the fit entirely results in similar 
central values for and /3jv, but with a larger statistical uncertainty. 

Finally, the constraint u, y\t\ < 160 MeV can be modified to uj, \f\t\ < 200 MeV, 
which allows all of the data to be fitted. The results of that fit are shown in Fig. an d 
are presented in Table El together with the three other fits discussed so far. 
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Figure 15: Results of the 0(Q 4 ) calculation at four different lab photon energies: 49, 
55, 67, and 94.5 MeV. Data as in Fig. I1UI The solid line is the 0(Q 4 ) result (including 
very-low-energy resummation) with oln = 11.5 x 10 -4 fm 3 and (3n = 0.3 x 10~ 4 fm 3 (Fit 
IV), while the dot-dashed line is the (parameter- free) 0(Q 3 ) calculation. 



The numbers differ from those given in our earlier publication jTHj for two reasons. 
First, there we mistakenly fitted the SAL data at a cm. energy of 95 MeV, whereas here 
we have done so at a lab energy of 94.5 MeV. This causes a relatively small change in the 
results for most quantities in Tabled but it does increase the central value for obtained 
with the Nijm93 wave function by 1.0 x 10~ 4 fm 3 over the value quoted in Ref. JHI- A 
larger effect occurs because there was an error in the program used to generate the 0(Q A ) 
results presented for 7d scattering in Ref. [TH]. The amplitude used there was missing 
a factor of i in its spin-dependent parts. This affects the 0(Q 3 ) results by less than 1%. 
However, at 0(Q 4 ) the omission of this factor modifies the interference between the A3- 
A 6 pieces of the single-nucleon amplitude and the 0(Q 4 ) two-body currents. Correcting 
this mistake increases the predicted cross sections by about 10-15% with respect to those 
published in Ref. [Oil, and those which appeared in the first version of this paper [51)] . 
We emphasize that this mistake was not made in the 7p calculation reported in Ref. [IB*] 
and discussed in Sees. [3] and |U above. 

As discussed above, cross sections depend on the choice of deuteron wave function 
at the 10-20% level. This change in the 7d cross section has a significant effect on the 
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Fit 


Wqvp 

VVOjVC 


VPTV-l OW-PT1 prcrv 






yN 


v 2 /d o f 




function 


resummation? 


below 


(10~ 4 fm 3 ) 


(10~ 4 fm 3 ) 




I 


NLO xPT 


No 


200 MeV 


13.6 


0.1 


2.36 


II 


NLO xPT 


No 


160 MeV 


15.4 


-2.3 


1.95 


III 


NLO xPT 


Yes 


160 MeV 


13.0 


-1.8 


1.33 


IV 


NLO xPT 


Yes 


200 MeV 


11.5 


0.3 


1.69 


V 


Nijm93 


Yes 


200 MeV 


16.9 


-2.7 


2.87 



Table 3: Results of our different 0(Q 4 ) EFT extractions of isoscalar nucleon polarizabil- 
ities from low-energy coherent 7d scattering data. 



extraction of and (3^. Changing only the wave function yields the difference between 
Fit IV and Fit V in Table El The high x 2 of Fit V is again due to a failure to reproduce 
the 49 MeV data — this time even when the very- low-energy resummation of Section l6~3*l is 
performed. While the Nijm93 wave function is not consistent with xPT, it does have the 
correct long-distance behavior. The differences we see are therefore a consequence of the 
different short-range behavior of the Nijm93 and NLO %PT NN potentials. As such they 
should be renormalized by 'jNN — > '-/NN contact operators. However, Weinberg power 
counting predicts that the appropriate operators do not appear until (at least) 0(Q 5 ). 
Thus, the degree of variability between the results with the NLO %PT and Nijm93 wave 
functions could be said to be inconsistent with Weinberg power counting. Indeed, it may 
be necessary to modify the Weinberg power counting so that such contact operators appear 
at a lower order than is indicated by naive dimensional analysis. Further understanding 
of this issue is clearly crucial to the use of %PT as an accurate calculational tool for low- 
energy reactions on deuterium, but lies beyond the scope of this paper. Here we adopt the 
conservative strategy of assigning a theoretical error which encompasses both the Nijm93 
and NLO %PT results. 

The x 2 contours for the five fits of Tableware shown in Fig. EE with the "best fit" EFT 
(Fit III) giving the 1-sigma region bounded by the solid black line. Putting these results 
together we conclude that our central values, and error bars, for the isoscalar nucleon 
polarizabilities are: 

a N = (13.0 ± 1.9)±f:| x 10~ 4 fm 3 , 

p N = (-1.8 ± 1.9)tH x 10~ 4 fm 3 . (53) 

The errors inside the brackets are statistical, and those outside reflect the arbitrariness 
as to which data are included, and which deuteron wave function is employed. 

The theoretical error should also include a component from higher-order (0(Q 5 ) and 
above) pieces of the 'yNN kernel. A number of these are taken into account in the 
sophisticated potential-model calculations of Refs. j2Hl 123 I2H1- Including them in the 
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Figure 16: Allowed regions (at the l-cr level) for isoscalar nucleon electric and magnetic 
polarizabilities (in 10 _4 fm 3 ) resulting from different 0(Q i ) EFT fits to low-energy data. 
In the order the fits are listed in Table H3 the constraints are: fit I (short-dashed, red); 
fit II (medium-dashed, orange); fit III (solid black); fit IV (long-dashed, green); fit V 
(dot-dashed, blue). For comparison, the Baldin sum-rule constraint found by averaging 
the proton and neutron results of Ref. (57] is shown as the gray shaded band. 

calculation will alter the central values, although the theoretical error arising from wave- 
function dependence and the variation in the kinematical cuts we impose is large enough 
that we expect it encompasses the impact of higher-order terms in the kernel. 

Combining Eqs. (jHOj) and (JH3)) we see that a wide range of neutron polarizabilities is 
consistent with a model-independent analysis of the current low-energy 7p and coherent 
7d data. However, there is clear statistical evidence in these data sets for a n being 
significantly larger than f3 n : (3 n is consistent with zero within our large error bars. Also, 
the results presented here provide no evidence for significant isovector components in a 
and P, in contrast to some previous results in the literature (see, e.g. Ref. [20J). 

6.5 Other techniques for measuring a n 

A dispersion sum rule — analogous to the Baldin sum rule for the proton — that relates 
the sum of isoscalar nucleon electric and magnetic polarizabilities to an integral of the 
deuteron photo-absorption cross section can be derived, if one assumes that the single- 
nucleon contribution dominates the deuteron photoabsorption cross section. Using this 
sum rule, and subtracting their result for a p + /3 P , the authors of Ref. jHZj found: 

a n + p n = (14.40 ± 0.66) x 10~ 4 fm 3 . (54) 

This result is not consistent with our "best fit" result if only statistical errors are con- 
sidered, but if one includes the theoretical errors in Eq. (|53[) there is no discrepancy. 
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Moreover, the theoretical error on the sum rule associated with assuming that the free 
neutron photoabsorption cross section can be found by taking the difference of deuteron 
and proton data is not clear. 

A more sophisticated evaluation of the neutron photoabsorption integral that uses the 
neutron multipoles of Ref. [SH] at low energy, together with experimental deuteron and 
proton photoabsorption cross sections at intermediate energies, and Regge phenomenology 
at energies above 1.5 GeV yields |27| : 

a n + p n = (15.2 ±0.5) x 10~ 4 fm 3 . (55) 

The theoretical error given in Ref. [2Zj is arrived at by employing different neutron pion 
photoproduction multipoles in the sum-rule integrand for the region up to ui = 500 MeV. 
However, the theoretical errors from the multipole analyses themselves seem not to have 
been assessed. Further, the extent to which the theoretical error in the high-energy part 
of the amplitude affects the final result is not clear to us. 

Most direct information on a n has been obtained by scattering neutrons on a heavy 
nucleus and examining the cross section as a function of energy. Controversy remains over 
what result this technique gives for a n . In Ref. [HE] the value 3 

a n = (12.6 ± 1.5 ±2.0) x 10~ 4 fm 3 (56) 

was obtained, which disagrees considerably with the result of the experiment of Ref. [U2] > 

a n = (0.6 ± 5.0) x 10~ 4 fm 3 . (57) 

Recent data from MAMI [21] exist for quasi-free 'yd — > ^pn in the 200 to 400 MeV 
range. These experiments update the pioneering work of Ref. [22], which was only able 
to set an upper bound on a n from this process. Ref. [21] reports the measurement of the 
reaction 'jd — *> 'jpn for laboratory photon scattering angles of 136.2 degrees, and then 
uses a theoretical model to extract: 

a n -fi n = (9.8 ± 3.6 (stat.))ti;l ( s Y st -) ± 2.2(model) x 10" 4 fm 3 , (58) 

from their data. This refines the constraint on a n — (3 n obtained somewhat earlier at SAL, 
where quasi-free Compton deuteron breakup was measured in the narrower range between 
236 and 260 MeV of photon energy and analyzed using the same theoretical model 
The result (p)8j) is consistent with the broad range for a n — (3 n indicated by our analysis 
of the low-energy coherent 7p and 7d data. 

3 The values measured in these experiments are not the polarizabilities (often denoted a) we have been 
discussing here. There is an additional term k„/M [HOI, with n n the magnetic moment of the neutron, in 
a n , which we have added to the experimental results from neutron-nucleus scattering in order to facilitate 
direct comparison. See also Ref. 
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7 Conclusion 



Starting from the most general Lagrangian that includes pions and nucleons and shares 
the global symmetries of QCD, we have calculated the amplitude for Compton scattering 
on the deuteron in xPT up to 0(Q 4 ). This amplitude consists of one- and two-nucleon 
pieces averaged over incoming and outgoing deuteron wave functions. Both the ~/NN 
kernel and the deuteron wave function were obtained from the same Lagrangian, with a 
consistent set of parameters. 

The one-nucleon piece of the jNN kernel is the same as the amplitude for Compton 
scattering on a nucleon. This was calculated to 0(Q A ) in the usual %PT power counting 
in Ref. That calculation contains four unknown parameters, which correspond to 

short-range contributions to electric and magnetic polarizabilities of the proton and the 
neutron. In this work we fitted the world's low-energy proton data f71Ei] and determined 
the two proton polarizabilities, obtaining values similar to the Particle Data Group's [oP] . 

We also calculated the two-nucleon kernel to the same order, using the power counting 
suggested by Weinberg [13J. It consists of various one-pion-exchange contributions with 
no unknown parameters. Finally, we used deuteron wave functions fully determined in 
the same approach [45J. To this order, the deuteron amplitude is sensitive only to the two 
isoscalar combinations of the nucleon electric and magnetic polarizabilities. From the 7d 
amplitude we calculated the differential cross section and compared it with deuteron data 
at various photon energies below 100 MeV ^IJJ 12H1 E] A good fit was obtained and we 
extracted from it results for the two isoscalar polarizabilities, which, together with our 
results for proton polarizabilities, allow an extraction of neutron polarizabilities. Various 
theoretical systematic uncertainties were studied. 

As pointed out in the Introduction, various sophisticated calculations of this process 
exist using potential models j2Hl 123 I2H1- Two of the distinctive features of our work are 
the use of a consistent xPT framework throughout, and the fact that we employ only low- 
energy data. We have used only data that we believe — on the basis of statistical tests — 
falls within the domain of the EFT. We find that a wide range of neutron polarizabilities is 
consistent with a model-independent analysis of this 7p and coherent 7d data. Narrower 
ranges for the neutron polarizabilities can be obtained from the data at the expense of 
additional assumptions which introduce model dependence. 

The potential-model calculations of Refs. 123 I2E] can be broken up into one- and 
two-body kernels and wave functions in the same way ours can. However, typically the 
different pieces are not calculated using a fully consistent field theory. The wave function 
is calculated from a phenomenological potential model that does not include, e.g., the two- 
pion exchange required by quantum field theory, but instead employs various single-meson 
exchanges. The best calculations strive for consistency of the 'yNN kernel with the two- 
nucleon potential by including in it the same mesons with the same form factors. This 
also leads to a gauge-invariant 7d amplitude if the necessary photon seagull diagrams 
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are included [27]. Meanwhile the 7 A" piece of the 7AA kernel in Refs. [25] [2TJ EE] 
includes a polynomical expansion of the 7N amplitude. They thus omit pieces of the 
7 A amplitude which are not analytic in the photon energy ^3]. On the other hand, the 
model calculations include a number of contributions that in our framework appear only 
at higher order. There is qualitative, but not quantiative, agreement between the different 
calculations and also between these potential-model calculations and the calculation 
reported on here. 

This is perhaps not surprising, since the potential-model calculations can be understood 
as various approximations to effective field theory. Indeed, it is one of the features of EFT 
that models which do not badly violate symmetries of QCD should — at low energies — 
reduce to the EFT with a particular choice of low-energy constants. (This choice may, or 
may not, agree with the choice that is made by QCD.) Contrary to what is apparently a 
common misconception in the literature, there are no more unknown parameters in our 
calculation than in these potential-model calculations. 

Various previous extractions of polarizabilities have also differed from ours in the em- 
phasis given to low-energy data. Frequently the Baldin sum rule [HI] — which relies on 
certain (reasonable) assumptions about the behavior of QCD amplitudes at higher ener- 
gies, and on model-based interpolation of data — is used. Our results are consistent with 
the Baldin sum-rule results for a p + (3 P [TKj . but do not rely on them. In addition, our 
results are consistent, within error bars, with the recent extraction of ± Pn from the 
Lund data using the detailed model of Levchuk and L'vov [211 ETJ- (But see also the 
values found using the data of Refs. 01201 and this model [2T>] l2"7|.). 

There are several ways to improve our present calculation. First, a complete treatment 
of the very-low-energy region would allow inclusion of low-energy data that could be 
measured at HI7S [55] . Second, one would like to systematically go to 0(Q 5 ) and higher 
in the present energy regime. The significant dependence of our results for a at and, to a 
lesser extent, on the choice of deuteron wave function also warrants further work. It 
could have important implications for the formulation of a fully-consistent nuclear effective 
theory. 

A promising future direction is the addition of an explicit Delta-isobar field, which 
would allow us to address higher-energy data using the same EFT philosophy employed 
here. A recent breakthrough in power counting has facilitated the use of ;^PT for reaction 
energies up to the Delta peak, in particular in the case of Compton scattering on the nu- 
cleon [HH] . Compton scattering on the deuteron has very recently been tackled along these 
lines, and it seems that the inclusion of an explicit Delta degree of freedom ameliorates 
the diffculties we had in reproducing the backward-angle SAL data [51] . 

Finally, there is, in principle, no obstacle to carrying out EFT calculations for the 
quasi-free process 7c? — > •ynp, and including the existing data [221 ESI 121] in the extraction 
of neutron polarizabilities. We hope that these various improvements will eventually lead 
to the use of EFT as a model-independent framework within which nucleon polarizabilities 
can be extracted from experimental data gathered in a number of different reactions. 
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A Third-order loop amplitude 



The third-order one-nucleon loop graphs give the following amplitude in a general frame. 
The notation tj is used for the tensor structures which multiply the amplitudes A± of 
Eq. (|2T)j): for example t\ = e' • e. 



T 
T- 

-L 7 7 



T 



VI 



ti + t 3 ) Jq[u, m^} + crossed 



vi 



2 1 rl 

9 A e ' 



V 



s 2 f 1 

— (ti + t 3 ) / dx( Jo[u - xuj, ml] + Jq[uj - xu', ml}) + crossed 

7T JO 

f ly r 

/ dx (d + l)ti Jq[ui, m\ — xyt] 
'o L 

+ (2(x + y — (d + 3)xy)ujuj't2 — 2V(x,y)ti — 2(\ — x — y)uu't^jQ\Cb,m 2 1T — xyt] 
+4u>u'xy(v(x, y)t 2 + (1— x— y)uouj't^j Jq[uj, m\ — xyt] 

+2(yia ■ (e" x q)e ■ k' + xia ■ (e x q)e' ■ k) J' q [6j, ml — xyt] 



dx 



(d — l)J Q [(x — l)Sco, ml — x(l — x)t] 



-2x{\ -x)(t- 5u 2 )J' [(x - \)5u, ml - x{\ - x)t] 



crossed 



(59) 



where u = (u + u)')/2, Cj = (UJ — xuj — you'), q = k — k' , t = (k 1 — k) 2 = 2uu'(cos9 — 1), 
5uj = uj — ijj'. The integrals Jo[u,rri 2 ), J2[u},m 2 ) and A^-fm 2 .] have their usual meanings P, 
prime denotes differentiation with respect to ml, and 

V(x, y)=u 2 + uuj' -ZJ 2 + \t{l -x-y + 2xy). (60) 
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We have also introduced an extra tensor structure: tj = a ■ {k' x k)t' -ke-k'. This is not 
independent; t 7 = sin 2 6*t 3 + cos #t 5 — 1 6 , but as it arises naturally in the calculations — and 
enters at too high an order in u to affect the polarizabilities — it is a useful notation. 

The notation "+ crossed" means that to every term is added another with e <-> e ', 
k <-» —k' and u —oj'. Since the tj are all either symmetric or antisymmetric under this 
transformation, the net effect is to add a term with u «-> —a;' to the coefficients of ti and 
t2, and subtract such a term from the coefficients of t 3 — t-j. 

In HB^PT the difference in photon energies, Sou, is of order Q 2 /M 4 . Thus for a 
fourth-order calculation only the expansion to first order in 5u> is required. This gives the 
following amplitudes 



,2 „2 



T 
T- 



Q 6 

-^-(ti + t 3 ) J [uJ, m 2 ] + crossed 



4^(ti + t 3 ) [ dxJ [. 

Jn JO 



xu, m J + crossed 



9_W 

VI 



dy / dx< (d + l)ti Jo [a;, m 2 — x?/t] 

n Jo 

+ (2(x + y-{d + 3)xyp 2 t 2 - 2V(x, y)ti 

-(x + y)cJ 2 (t 6 - t 5 ) - 2(l-x-y)uJ 2 U^j Jop,m 2 - xyt] 

+§(x - y) 2 5uj( id ■ (e" x q)e- k' - iff ■ (ex q)e' • 9J ^ J m " ~ ^ 



duj 



+Axyw 2 (V(x, y)t 2 + (1 — x — y)u) 2 t 7 ) Jq [cu, m 2 — xyt] j + crossed 



T 



VI 



dx 



(d- l)J [0,m 2 -x(l -x)t] 

-2x(l - x) t Jq[0, ml - x(l - x)t] 



where a; = (1 — x — y)u, and 



V(x, y) = oj 2 + ht(l — x — y + 2xy). 



(61) 



(62) 



There is only one term proportional to 5uj (in Ta 
the fourth-order amplitude (see below). 



and this cancels against a piece from 



4 Strictly speaking, all of the energy arguments of the Jo's should include v -p+, with p + the average of 
the initial and final heavy-baryon (off-shell) four-momentum. This effect, arising from the nucleon kinetic 
energies, alters the amplitude only at order Q 4 , and has been suppressed in the expressions above. The 
contribution is included in the fourth-order amplitudes calculated in Ref. and used here. 
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B Fourth-order loop amplitude 

The full amplitude in the Breit frame for diagrams |2M2fc can be obtained from Ref. I n 
a general frame, there are, in addition to the Breit frame terms, contributions containing 
p + , the average of the incoming and outgoing nucleon momenta. It is useful to define the 
following tensor structures 

t' 2 = e' ■ p + e ■ k' + e ■ p + e' ■ k 

t' 4 = id ■ (p + x q) e' ■ e 

t' 5 = id ■ (e' x k) p + ■ e — id ■ (ex k') p + ■ e" 

tg = ia ■ (e' x k') p + ■ e — id • (ex k) p + ■ e" 

tg 6 = id ■ (e' x p + ) k' ■ e — id ■ (ex p + ) k ■ e' 

t' 7 = id ■ (k' x k) (e" • p + e ■ k' + e ■ p + e' ■ k) 

t' 7 ' = id-{p+xq)e' -ke-k' (63) 

The notation reflects the structures that these t's give rise to in the center-of-mass frame. 
Again, not all the structures are independent, but it greatly simplifies matters to work 
with this set. 

Diagram by diagram, the p + -dependent contributions are 

a = ~2MP^ l+ 3)P+ + Thj + crossed 

T b = ^(2(t 1 + t 3 )p + -fc + + (t / 2 + t / 5 )) fdxx dM g2 ml] + CT ossed 

T f = — ~ ^ (l - r 3 ) t' 6 c^ 1 / cfeJofxu;, m%] + crossed 
4iW Jtt Jo 

T g = ~ 4 ^ 2 (1 + r 3 )t 6 u;~ i y rfxJ [x^,m^] + crossed 



2 2 /-l rl-y 
Th = -Mt^ dy dx 



2Mf, 
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v l-ar-2/)p+ • £;+ ( ti(d + 3) — 



- ((x + yV 2 (t 6 - t 5 ) + 2{l-x-y)u% + 2V(x, y)t a 

\ dJ' [Co,ml - xyt] 



+Axyu 2 (V(x, y)t 2 + (l-x-y)u 2 t 



+2((d + 5):n/ - sc - 2/)ort 2 J 
dJo[u,ml - xyt] 



duj 



- (|((d + 3)(x + y) - 2)t' 2 - < + t' 6 - tf 5 ) 

+ ((1 -x-y) (x + y)co% - 2xycuH>; + (x + y) y)tQ g^gi^Z^! 
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1 -5a; / / dx{x - y) i(a ■ I x qe ■ k' - a ■ e x qe' ■ k) 1 1 
Jo Jo 



+crossed 

r< = -f^tf 4 fA;K-x(i-x)t] 

2 g 2 /-l 

T i = 177^ / rf W rfx 2^t^ 2 A';K - a^t] - ~ ^ + crossed 

M /-7r JO JO L 

T k = -Ti (64) 

Here we have defined k + as the average of the incoming and outgoing photon momenta. 
The tensors t[ have the opposite crossing symmetry to the corresponding tj, as do the 
products p + ■ k + ti. 

Note that the piece in proportional to 5uj comes from p + ■ q, and exactly cancels the 
piece proportional to 5uj in the expansion to fourth order of the third-order amplitudes. 

All the rest of the terms written above can alternatively be generated by a boost of 
the third-order Breit-frame amplitude. The transformations to do this are as follows: 

with analogous transformations for the other photon. These boosts are odd under crossing, 
which accounts for the change in symmetry mentioned above. It is easily seen that the 
pieces of (J64J) proportional to p + ■ k + come from the boost of u; they are generated from 
the fourth-order diagrams which are just third-order diagrams with an insertion from 
on the nucleon line. The other structures come from the boost of the vectors e, k, etc. 
Applying the boost of Eq. (J55)l in the third-order amplitude gives the additional terms 
needed to ensure equality with the result from explicit evaluation of the fourth-order 
diagrams in an arbitrary frame. 



C Born terms 

In addition to these terms arising from pion loops in a boosted frame, there are also Born 
terms. The complete set of Born terms corresponds exactly to the expansion in powers of 
1/M of the result obtained from Dirac nucleons, with the vertex of an incoming photon 
of four-momentum q being 

r M = Z llx + —a^q", (66) 
as long as loop renormalizations of bare parameters such as k are ignored. 
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Born terms fall into two categories. First, there are those which do not vanish in the 
Breit frame, 

2 

TiZt = —^\((Z + K) 2 (l + cos9)-Z 2 )(cos9-l)tj 2 t 1 -K(2Z + K)cos9ujH 2 
4M d LV / 

2 ^2 

' ~ ;E^t x (67) 



2M 2 

where is the average kinetic energy of the nucleons. The coefficient of the nucleon 
kinetic energy term depends on the normalization of the Dirac spinors in the relativistic 
case. As shown, the normalization is uu = 2E/(E + M), which is appropriate for the 
standard non-relativistic reduction. However, the usual covariant treatment of the Comp- 
ton Born terms corresponds to the normalization uu — 1, in which case the last term on 
the right-hand side of ()67|) is not present. 

Second, in an arbitrary frame, there are the terms generated by a boost of the LET 
pieces of the third-order, Breit-frame amplitude 

e 2 

T-iBorn ° -?2 + / 

boost - M2u * h 



e 2 



p + ■ k + [{Z + /0 2 (t 3 (cos# - 1) + 1 4 - t 5 ) - /t 2 t 3 + Z{Z + «;)t 6 ) 



2M 3 

+ (^Z 2 - [Z + k) 2 ) t' 4 + k% - Z(Z + K) (1 - cos 6% + k{Z + K )t' 56 
+ ^Z 2 p + ■ k + t' 2 + 2Z 2 (cos 9 - l)e • p + e ' ■ p+ . (68) 

The first part of the t' 4 term comes from the Thomson term via a Wigner rotation, which 
is a second-order boost effect; the terms in the last line are also second-order, with the 
second-order boost given by 

e-p + ~ (p + ■ k)(e-p+) £ e-p + _ 

k - % -^V+ + P ~i^P+ ( 69 ) 

The angle dependence of these Born terms is more complicated than that associated with 
the polarizabilities. 

D 7^V kinematics in the 7c! center-of-mass frame 

In this appendix we consider the evaluation of the single-nucleon Compton amplitude 
"inside" the deuteron. This necessitates the evaluation of expectation values: 

wm, (70) 
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where the t[ are the structures defined by Eq. (JHHJ), and \ip) is a deuteron wave function, 
calculated for a deuteron at rest. In the next appendix we will explain how to correct 
such matrix elements for the fact that the initial-state (final-state) deuteron wave function 
must be calculated in a reference frame which is moving with momentum —k (—&'). 

All t^s depend linearly on p + , the average of initial- and final-state nucleon momentum 
in the frame of choice. With the kinematics defined in Fig. [7| this average, in the frame 
where the jN collision takes place, is: 

p + = l(p + p')-i(k + k'), (71) 

with p' = p + i<f . 

The evaluation of the expectation values ()70|) then involves integrals of the form: 

(1 " p r(p')pip(p), (72) 



(27T) 



where all spin labels are suppressed. But, for deuterium, which contains only L = and 
L = 2 components ip(p) — ip(—p), which implies that: 

<rp rip') (p + p')4>(p) = 0. (73) 



(27T) 

Combining (J75|) and (J7TJ) we see that the expectation values (fTUj) may all be evaluated by 
making the replacement: 

p + -,-\(k> + k) (74) 

This substitution is, however, only valid because we are considering 7d elastic scattering. 

Now, for the third-order 7d amplitude none of the boosts of the third-order loops are 
needed. The only effect arises from boosting the Thomson term, using the expressions 
(jfi5j) for the boosts of the polarization vectors. (An equivalent result is obtained from 
computing the u-channel nucleon pole with El ^yNN vertices.) This yields the piece 
of the 7iV amplitude which appears on the first line of Eq. (|68|h and so the only t[ 
whose expectation value is needed at 0(Q 3 ) is t' 2 . In the 7d center-of-mass frame this 
contribution to the 0(Q 3 ) amplitude can be evaluated using the replacement (|74*jl . 
The sole effect of such a boost is a modification of the function A2. 

Af c - m ' = A* rcit + (75) 
2 2 2M 2 v ; 

Explicit numerical evaluation of the full expectation value of t' 2 shows that the replacement 
(|74|) is accurate to a very good approximation. Corrections to it arise from terms in the 
deuteron wave function suppressed by uj 2 /Mj. Nominally these are 0(Q 5 ), but in reality 
they are smaller still, since they have Mj in the denominator. 
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At 0(Q 4 ) we must also consider expectation values (fTUj) for i = 4,5, 6, 56, 7, as well as 
the expectation value of t" and p + ■ k + . The pieces of the general-frame amplitude which 
arise from boosts of the third-order amplitude were listed in Appendix [B] and on the 
second and third lines of Eq. (jBEJ). They are all linear in p + . Up to corrections of an order 
higher than we consider in this work, their contribution to (J7U|) in the 7d center-of-mass 
frame may be obtained via the replacement (J74)) . This produces the following results for 
A] d c - m \ i = 1...6: 



a ■yd cm. 


^Brcit 


+ 


CO 

2M 


j^jd, cm. 


^Breit 




uj 

2M 


j^yd cm. 


^Breit 


+ 


CO 

2M 


a ■yd cm. 
A 5 


^Breit 


+ 


CO 

2M 



[(1 -cos#)Af rcit -Af eit ] ; 
sin 2 9(Af eit + Af reit ); 



Af cit + 



e 2 co\ 
8M 3 ' 



1 - - cos 6 ) Af clt - cos 6{Af Clt + Al 



Breit 
3 



.4 



Breit \ 
6 J 



A 



7d cm. 
6 



a Breit , 
Si 6 + 



CO 
2M 



1 

2 



^Breit + | j _ - cog j ^Breit + - y4 Brei 1 



3 

2' 



(76) 



The boost of the Thomson term must be taken care of separately, since here the second- 
order boosts of Eq. (jQQj) are required. This leads to the fourth-order contributions to the 
7iV amplitude listed on the last line of Eq. (|68|). For one of these terms the expectation 
value 

(tp\e-p+ e'-p+\ip) (77) 

enters. It cannot be calculated by symmetry arguments, but must be retained and com- 
puted by explicit numerical integration. These 0(Q A ) terms are included in our calculation 
of the 7d amplitude, although since they are suppressed by 1/M 2 relative to leading their 
contribution is numerically small. 

Lastly, in the general-frame, 0(Q A ), Compton amplitude discussed in AppendiceslHland 
[HI there are a number of terms proportional to p + -k + . These terms all arise from applying 
the first-order boost (|63|) to the Breit-frame photon energy ujf,. Using the replacement 
()74j) — which is valid for the terms linear in p + — we may replace p + ■ k + by — \k + \ 2 /2. 
Consequently all but one of the these terms in the fourth-order amplitude may be included 
in the final result simply by employing the energy 



1 



^halfb = 



AM 



CO (1 + cos( 



(75 



when evaluating the third-order 7iV amplitude. This is the strategy we adopted in our 
calculation. 
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. only term in the 0(Q A ), general- frame, 'jN amplitude for which the use of coWfb m 
the 0(Q 3 ) amplitude does not yield the correct p + ■ k + pieces of the 0(Q 4 ) amplitude 



The 

3 0((, , 

the penultimate term in Eq. (|68|): 

- ■ ^ (79) 

which is second-order in p + . The error that results from evaluating this piece of 
using u ha _i{b in the third-order amplitude is of relative order 1/M 2 . Ultimately, we expect 
a number of 0(Q 5 ) effects to be more important than this particular difference. 

E Boosting the deuteron wave function 

At 0(Q 4 ) we also have to take into account that the deuteron is not in its rest frame in 
either the initial or final state. 

Consider the two-body matrix element of the operator O: 

(0} cm = (p', -k'\6\p, -k), (80) 

where we are employing a basis of two-nucleon states expressed in terms of relative and 
center-of-mass momenta, and we have chosen to work on the 7d center-of-mass frame, 
where the initial (final) momentum of the two-nucleon system is — k (—k'). Using the free 
boost operator xo and working to first order in the boost (which is all that is necessary 
at this order) we write: 

(0) cm = (0) rest + (p', O\i\xo(-k')0 - dxo(-k)]\p, 0). (81) 

where: 

(d)rest = (f, 0\O\p, 0) (82) 

is the zeroth-order result. 

The free boost operator Xo(P) is given by 



X( p ) = [~\ [ f - Pp-P + p-Pr-P]+^i- e? 2 ) xp-Pj, (83) 

where r and p are to be interpreted as quantum-mechanical operators, i.e. they do not 
commute with each other. So we see that we can write 

(0) cm = (d) rest + (d) jCr + (d) ±a , (84) 

where, upon evaluation: 

(6)xr = ^{ U}2+ P > - V p' -k' + p-kV p -k) (6) rest , (85) 
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and, if O is spin-independent: 



(6)x« = g^2 [i S = ~ *») ■ fi' xp'-kx p)\S = 1>] (6) rest , (86) 

which vanishes because the \S = 1) wave function is symmetric under the interchange of 
the spins of particles one and two. 

Since both {0)% T and (0)% a are suppressed by 1/M 2 at 0{Q A ) we need only consider 
what impact they have on the evaluation of the 0(Q 2 ) ^NN kernel. Thus, we now move 
to the specific case where the operator O is the 0{Q 2 ) 7A amplitude: 

(6) rest = -^(p'-p-\{k-k>)). (87) 
After some algebra we find that: 

(O)rest + (0), r = (l + ^) ^ (P'-P- Weff ~ Kf f)) , (88) 

where: 

jg e// = jg + g -^ + /) fc-A4fc, (89) 

n 8M 2 16M 2 ' v ; 

with a similar result for k' e ft- Using the fact that p' is constrained to be p + \q when we 
are evaluating this particular operator we find: 

ft// " *=// = * - wfi - » s «>) ~~ k+k IKF * ~ U l < 90 » 

If we write p as — plus terms which vanish upon taking the expectation value at the 
order we work to here, Eq. becomes: 

^//~^/ = ? ( 1_ T^) ' (91) 

which coincides with Adam and Arenhovel's result for the one-body charge operator in 
the Breit frame jHS]. This effective reduction in k' — k can be interpreted as a length- 
contraction effect. However, at the energies considered here it is very small: at 9 = 180 
degrees, uj = 95 MeV, it produces only a 0.3% change in \q\. 
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